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FOREWORD 


This report was prepared b/ Aerospace Systems, Inc. (ASI), Burlington, 
Massachusetts, for the National Aeronautics and Space Administration (NASA) under 
Contract No. NAS 2-8671 , The report documents the results of research performed 
during the period October 1975 to May 1977. The study was sponsored by the Ames 
Research Center, Moffett Field, California, with Dr. Wayne Johnson serving as 
Technical Monitor. 

The effort was directed by Mr. John Zvara, President and Technical Director 
of ASI. Mr. Richard B. Noll, Vice President of ASI, served as Project Engineer. 

Mr. Paul Soohoo of ASI and Dr. Luigi Morino, Director of the Computational Continuum- 
Mechanics Program at Boston University, were Principal Investigators. Dr. Norman D. 
Ham, Director of the V/STOL Technology Laboratory at MIT, contributed to the study 
as Principal Consultant, 
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SUMMARY 


An investigation has been performed to demonstrate the feasibility of the 
Green's function method to study tilting proprotor aircraft aerodynamics with particular 
application to the problem of the mutual Interference of the wing-fuselage-tail-rotor 
wake configuration. While the formulation is valid for fully unsteady rotor aerodynamics, 
this report is mainly concerned with steady-state aerodynamics, which is achieved by 
replacing the rotor with the actuator disk approximation. The use of an actuator disk 
analysis introduced a mathematical singularity into the formulation; this problem has 
been studied and resolved, 

Tbe pressure distribution, lift, and pitching moment were obtained for an 
XV-15 wing -fuse I age -tail rotor configuration at various flight conditions. For the 
flight configurations explored in this investigation, the effects of the rotor wake inter- 
ference on the XV-15 tilt rotor aircraft yielded a reduction in the total lift and an 
increase in the nose-down pitching moment. 

The present method has potential application in the design of helicopters 
because it provides an analytical capability that is simple to apply and can be used to 
investigate fuselage-tail rotor wake interference as well as to explore other rotor design 
problem areas. The numerical results obtained are of a preliminary nature and were 
used to identify requirements for further theoretical and experimental research into 
wing -tail -rotor aerodynamic interference effects on tilting proprotor aircraft. 
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SECTION 1 


INTRODUCTION 

The XV-15 tilHng proprotor (tilt rotor) aircraft is being developed under a 
joint NASA/Army program to provide flight research and technical data for the design 
of tilt rotor aircraft for military and civil applications (References 1 and 2). This tilt- 
ing proprotor research aircraft is characterized by low disk-loading rotors (two 3-bladed, 
25-ft diameter rotors for a gross weight of 13,000 lb) mounted at the wing tips, a forward- 
swept high wing, and twin vertical tails. The XV-15 aircraft offers the hover and low- 
speed capabilities of a helicopter and the high-speed cruise (300 knots) advantages of 
fixed-wing aircraft. The proprotors are positioned in a helicopter mode to take off 
vertically. As forward speed is developed, the rotors are tilted forward progressively 
(conversion mode) until the high speed cruise mode is achieved. In this aircraft mode, 
lift is provided by the wings, and the proprotors provide the propulsive thrust. 

The large rotors located on the wing tips and tilted at different thrust angles 
introduce numerous dynamic considerations (e.g.. References 3 to 6). In particular, 
the aerodynamic wake generated by the rotors can envelope the tail surface. Thus, 
the mutual interference of the rotor, wing, fuselage, tail, and their wakes *s poten- 
tially of greater Importance for this aircraft than for other configurations, (compared to 
helicopters, the aircraft has larger wing and tail surfaces; compared to propeller air- 
planes, it has much larger rotors which perform the same function as propellers. The 
effects of the aerodynamic interference on the stability and control characteristics of 
the aircraft are expected to be important. Of particular interest is the effect on longi- 
tudinal trim of the aircraft, especially through transition and conversion modes. 

This report demonstrates the feasibility of a method to study aerodynamic 
interference effects on aircraft. In particular, the problem of the rotor wake interaction 
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effects for a tilting proprotor aircraft at various flight conditions is investigated. The 
results obtained are of a preliminary nature and are used to identify requirements for 
further theoretical and experimental research into wing-tail -rotor aerodynamic inter- 
ference effects on tilting proprotor aircraft. 

The report begins in Section 2 with o description of the tilt rotor aircraft and 
a review of problems anticipated because of aerodynamic interference. The state of the 
art in the prediction and Oiialysis of aerodynamic interference is reviewed briefly in 
Section 3. The approach to the problem used in this study is discussed in Section 4, along 
with required modifications to a computer program used in calculating aerodynamic loads. 
Section 5 presents the potential aerodynamic theory for the mutual Interference of a wing- 
fuselage-tail-rotor wake configuration. In Section 6, the potential aerodynamic theory 
is modified for the special cese of an actuator disk model ! istead of o finite-blade rotor 
model. Section 7 gives the numerical formulation used to implement the aerodynamic 
theory of Sections 5 and 6 into a computer program. The determination of appropriate 
flight conditions to be used in the investigation is discussed in Section 8. Section 9 
presents the computed results and discusses the implications of the calculations. Section 
10 provides concluding remarks on the present study, and gives recommendations for 
future research based on the results of this study. A list of references supporting the 
various discussions is included, as well as appendices which provide supplementary 
information on the computer program used during the study and on the actuator disk 
equation. 
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SECTION 2 

TILTING PROPROTOR AIRCRAFT BACKGROUND 
2.1 AIRCRAFT DESCRIPTION 

The tilting proproior aircraft combines in one aircraft the hover efficiency 
of ;'ie helicopter and the high-speed efficiency of the fixed-v/ing aircraft (References 
1 7, and 8). The NASA/Army tilting proprotor research aircraft program sponsored 

pniliminary design studies (References 9 and 10), which led to the development of a 
Hit-rotor research aircraft designated as the XV-15. 

The XV-15 (Figure 1 and Table 1) is representative of those aircraft which 
en )loy the tilt rotor concept. The hover lift and cruise propulsive force is provided by 
iovif disk -loading rotors located at each wing tip. The rotor axes rotate from the verti- 
cal (or near vertical), the normal position for hover and helicof>ter flight, to the 
horizontal for air'Iane mode flight. Conversion characteristic: are given in Table 2. 
Hover control is provided by rotor generated forces and moments, while airplane mode 
flight control is obtained primarily by the use of conventional aerodynamic control 
surfaces (Table 3). 

A c ss-shafting system connecting the rotors provides several benefits. 

This system precludes the complete loss of power to either rotor due to a single engine 
failure, permits power transfer for transient conditions, and provides rotational speed 
synclironization. Rotor-axis tilt synchronization is ccl.ieved by a separate interconnect 
shaft. 


aircraft is capable of high duration hover (approximately one hour at 
f ' design gross weight), helicopter mode flight, versatility in performing conversion 

I (steady-state flight is possible at any point in the broad transition corridor), and air- 

? 

i plane mode level flight at speeds greater than 300 knots. The low disk-loading rotors 

f. 

f 
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Table 1. Dimensions and General Data for XV-15. 


Description 

Wing 

Horizontal 

Tail 

Vertical 

Tail 

Area, sq ft 

169.0° 

50.25 

50.5 

Span, ft 

32.17° 

12.83 

7.68 

Chord, ft 

5.25 

3.92 

4.09/2.40^ 

c, ft° 

5.25 

3.92 

3.72*^ 

Aspect ratio 

6.12 

3.27 

2.33 

Incidence, degree 

3.0 

0 to +6® 

— 

Dihedral, degree 

2.0 

0 

— 

Sweep y , degree 
4 

-6.5 

0 

31.6^ 

Section (NACA) 

64A223 

(Mod) 

64A015 

0009 

Tail length, ft 

— 

22.4 

23. 2^^ 


a - between rotor centerlines at lilt axis; b - root/tip; c - mean aero- 
dynamic chord (MAC); d - total; e - hinged at 33 percent chord; 
f - upper. 


can be operated in the autorotation state to reduce descent rote in the event of total 
power loss. Research operation at the design gross v/cight allows for a total useful load 
of over 2900 pounds. At intermediate rotor-axis tilt angles (between 60° and 75°), the 
aircraft can perform STOL operations at weights above the maximum VTOL gross weight. 


The Lycoming LTC1K-4K (a modification of the r53-L-13B) engines and main 
tronsmissions are located in wing-tip nacelles to minimize the operational loads on the 
cross-shaft system and, with the rotors, tilt as a single unit. The use of the free turbine 
engines permits the reduction of rotor rotational speed for airplane mode flight to improve 

rotor performance and reduce cruise noise. 
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Table 2. Conversion Characteristics for XV-15. 


Conversion axis: 

Wing chord location, percent 

39.0 

Forward sweep, degree 

5.5 

Dihedral, degree 

3.0 

Angle-to-mast axis, degree 

95.5 

Angle of outboard tilt-of-mast axis 

Helicopter mode, degree 

2.5 

Airplane mode, degree 

0 

Conversion range, degree 

0 to 95 

Rotor tip clearance (helicopter mode) to: 

Ground, ft 

11.5 

Wing upper surface, ft 

4.8 

Rotor tip clearance (airplane mode) to: 

Fuselage, ft 

1.0 

Wing leading edge, ft 

0.47 


Table 3. Dimensions of Movable Aerodynamic Surfaces for XV-15. 


Description 

Flap 

Flaperon 

Elevator 

Rudder 

Area aft hinge, sq ft 

11.0° 

20.2° 

13.0 

7.5° 

Span, ft 

4.25*^ 

7.86^ 

11.0 

4.66 

Chord aft hinge, percent 

0.25 

0.25 

0.30 

0.25 

Surface deflection, degree 

c 

c 

±20 

±20 

Control travel for total 
surface deflection, in 

— 

9.6 

9.6 

5.0 


a - total both panels; b - one side, along hinge; c - see Figure 3. 
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The XV-15 research aircraft utilizes 25-ft-diameter gimbal-mounted, 
stiff-inplane, three-bladed rotors, with elastomeric flapping restraints for increasing 
helicopter mode control power and damping. Rotor^lade characteristics are shown 
in Figure 2 and Table 4. Variable rotor-tip speed control is provided to enable 
research on noise, performance, and hover downwash. The nominal design tip speeds 
are tabulated as follows: 


Condition 

Tip Speed, ft/sec 

rpm 

Hover/Helicopter Mode 

740 

565 

Cruise/Airplane Mode 

600 

458 

Hover Test Overspeed 

818 

625 

Design Limit 

865 



661 


The forward - swept wings provide clearance for 12° of flap deflection which will accom- 
modate gusts and maneuver excursions while operating in the airplane mode. Wing- 
pylon-rotor stability is accomplished by selei-ting a stiff wing and pylon-to-wing 
attachment and minimizing the distanc? of the rotor hub from the wing. Airplane 
mode wIng-pylon-rotor stability is retairied up to airspeeds of 370 knots even with a 
20 percent reduction in wing and pylon stiffness. 

For hover flight, the wing flaps and flaperons are deflected downward to 
reduce the ving download and incrt-se hovering efficiency (Figure 3). Hover roll 
control is provided by differential rotor collective pitch, pitch control by cyclic 
pitch, and yaw control by differential cyclic pitch. Pilot controls in the helicopter 
mode are similar to that of a conventional helicopter. A collective stick provides 
power and collective pitch for height control and a control stick provides longitudinal 
and lateral control. 
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Number of blades per rotor 
Diameter, ft 

Disk area per rotor, sq ft 
Blade chord (see Figure 2), in 
Blade characteristics 
Solidity 

Hub precone angle, degree 

pitch flap coupling, degree 

Flapping design clearance, degree 
Blade Lock number 


25.0 
491 

14.0 

See Figure 2 
0.089 
2.5 
-15.0 
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Figure 3 Flaperon Deflection versus Flap Position. 


In the airplane mode, conventional airplane stick and rudder pedals are 
employed, while the collective stick/power lever continues to be used for power 
management. An H-tail configuration (two vertical stabilizers) has been selected to 
provide improved airplane directional stability a»'5und a zero yaw angle. Control 
authority for the power lever, blade pitch governor, cyclic, differential cyclic, 
differential collective, and flap/flaperon relationship are phased with mast angle by 
mechanical mixing linkages. 


2.2 ANTICIPATED PROBLEMS DUE TO AERODYNAMIC INTERFERENCE 

The tilting proprotor is a VTOL aircraft with laterally displaced tilting 
rotors at the tips of small wings, and a conventional tailplano mounted at the rear of 
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the fuselage for stability and control in transintional flight (see Figure 1). The mutual 
interference of the rotor, wings, fuselage, and tail is likely to have important effects 
on the stability and control characteristics of this aircraft, including the longitudinal 
trim during transition and conversion, and the lateral trim in a sideslip or crosswind. 

In a previous aircraft of this type, the Bell XV -3, a nose-down pitching 
tendency was caused by the rotor woke inboard tip vortices passing near the horizontal 
stabilizer. At low airspeeds, the tip vortices passed below the horizontal stabilizer and 
induced an increase in the angle of attack. The resulting lift required a small aft 
stick displacement which reduced the stick displacement gradient in the forward speed 
range from hover to 40 knots. 

Tilt rotor wake effects on the wing and horizontal stabilizer have been mea- 
sured during a recent one-tenth scale powered-model wind-tunnel test (Reference 11). 
From this test it was possible to establish these effects during low-speed transition and 
conversion flight. Wind-turnel data showing the percent wing lift as a function of 
airspeed and mast angle are given in Figure 4. Comparison with rotors-off tests indi- 
cated that above 40 knots, v/ing lift is not influenced by the rotor wake. 

The tilt rotor wake effects on the horizontal stabilizer were determined 
from the data summarized in Figure 5. As shown, an upwash exists during low-speed 
flight. This was found to have a considerable effect on longitudinal trim. 
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Figure 4. Rotor/Airframe Interference in Forward Flight (Reference 11). 
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Figure 5. Rotor/Horizontal Stabilizer Wake Effects in 
Forward Flight (Reference 11). 
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ScCTION 3 
STATE OF THE ART 
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3.1 WING-BODY AERODYNAMIC INTERFERENCE 

The classical approach for evaluating the aerodynamic pressure on interfering 
wings and tails Is the lifting surface method (References 13 and 14). Calculations have 
been normally limited to wings or other lifting surfaces separately, except for a limited 
number of cases in which the vehicle could be considered as a slender body. However, 
interaction of the flowfields of v/ings and the fuselage results in interference effects 
(■Reference 15) v ’ ich cannot be adequately accounted for when examining the lifting 
surface and the body separately. The advent of modern high-speed computers has 
prompted the development of a series of new methods, generally called finite element 
methc is. These methods offer both flexibility and efficiency and can be applied to 
complex configurations. AH the methods are based on potential aerodynamics and 
consider various distributions of sources, doublets, vortices, pressure panels, etc., 
on and/or inside the surface of the aircraft. 

An early work on the flow field around three-dimensional bodies (Reference 
16) uses constant strength source -elements to solve the problem of steady subsonic flow 
around nonlifting bodies. This method has been extended to lifting bodies (References 
17 through 19) by including doublet and vortex panels. Woodward's method (Reference 
20) is a different approach which uses lifting surface elements for the representation of 
the body. For oscillatory subsonic aerodynamics, extensions of the doublet-lattice 
method (Reference 21) are obtained by either placing unsteady lifting surface elements on 
the surfoce of the body or by using the meihod of Images combined with slender body 
theory (Reference 22). In the supersonic unge, complex configurations are analyzedin 
References 20, 23 and 24. Recently, Woodward's method (Reference 20) has been 
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extended to provide u unified method for both subsonic and supersonic steady flow 
(Reference 25). In the wing -body -to 1 1 interference methods mentioned above, none 
is sufficiently general to analyze the interference of complex aircraft configurations 
in the presence of a rotor, since there exists no frame of reference with respect to 
which all the surfaces (body and rotor) are fixed. 

A general method for steady and unsteady subsrmic and supersonic aero- 
dynamics around complex configurations (where the surface is allowed to move with 
respect to the frame of reference) has recently been developed by Dr. Luigi Morino. 

Morino's technique (References 26 and 27) is based upon the Green function method 

for the equation of the aerodynamic potential. The potential is expressed in terms of 

source distributions with known intensity (obtained from the boun Jary conoitions) and 

doublet distributions with intensity equal to the unknown value of the potential on 

the surface. The integral equation is obtained by imposing the continuity of the 

potential along the normal to the surface; it may be noted that the boundary conditions , 

are automatically satisfied. The application of Morino's method to finite-thickness 

i 

wings is given in Reference 28, and the general theory for subsonic flow is reviewed 
in Reference 29. This method is readily applied to supersonic flow and is the only 
known method for analyzing unsteady supersonic flow around non-zero thickness con- 
figurations (Reference 3C). 

3.2 TILT ROTOR AIRCRAFT AERODYNAMIC INTERFERENCE 

There is relatively little information available on aerodynamic interference 
effects in tilt rotor aircraft. References 31 , 32, and 33 describe the results of tests of a 
powered wind tunnel model. Aerodynamic interference in forward flight was investigated 
by the observation of flow in the vicinity of the wing and empennage and the measurement 


4 

/ 
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of the effect of interference on the empennage. Smoke Injected upstream of the left- 
hand rotor showed the nature of the rotor-wing aerodynamic interference. Tests of a 
l/lO-scale model had previously shown that the wing download In hover changes to an 
upload in forward flight. The wing begins to contribute lift at cn airspeed of 35 knots, 
whereas linearized wake theory would predict the wing to be immersed in the roi .1 
wake at that low an airspeed. Observation of the flow showed that at speeds as low 
as 20 knots the rotor induced a strong upwash at the wing leading edge and that the 
rotor '/rake was completely off the wing at speeds over 30 knots. The nonuniform 
induced velocity of the rotor appears responsible for the wake moving off the wing at 
low airspeeds. 

As a result of the powered model tests, the following were determined: 

o Rotor wake interaction on the wing is nonlinear, varying from 
downwash on the outboard panels under the rotor to upwash on 
the inboard panels. 

« Rotor wake on the horizontal stabilizer varies from an upwash 
during low-speed helicopter flight to a downwash during high- 
speed flight. 

• The dynamic pressure at the empennage increases significantly 
during lo’^-speed flight to nearly double that of free stream. 

• The rotor wake across the horizontal stabilizer is nonlinear with 
sideslip and causes a pitch up with sideslip. 

• Directional stability is influenced by the rotor wake interoction or, 
the vertical stabilizer. 

The flow patterns in the vicinity of the empennage were studied with the 
aid of smoke and photographs of a tufted grid. At 30 and 40 knots, the flow patterns 
became distinct and showed that the rolled-up rotor wake vortices were located 
directly above the horizontal stabilizer. At these speeds, the rotor wake has the 
characteristics of the wake of a low aspect ratio wing, with the net effect of an 
upload at the horizontal stabilizer. 
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The tufted grid was not used at airspeeds higher than 40 knots. Smoke 
patterns at higher speeds showed that the rotor vortices moved downward and outward 
as speed increosed. A net upwash over the ^an of the horizontal stabilizer was visible 
at airspeeds between 60 and 100 knots. 

When the model was yawed at low speeds, the vortex cores shifted laterally 
with respect to the fuselage centerline. This shift has two effects on the empennage: 

(1) it reduces the strength of the rotor-wake-induced lift on the horizontal stabilizer, 
reducing the upload and producing a nose-up pitching moment; and (2) for yaw angles 
between 0 and 12 degrees, it reduces directional stability, due to the vortices moving 
laterally with respect to the vertical fins. At yaw angles greater than 12 degrees, the 
rotor wake increased the effectiveness of the vertical fins and thereby increased the 
directional stability over that of the basic airframe. 

p The interference between the rotor wake and the empennage was extracted 

>■ from test data under the assumption that the aerodynamic interference is superimposed 

j on the aerodynamic characteristics of the basic airframe. The downwash at the hori- 

zontal stabilizer was assumed to be caused by the wing and the rotors, with each acting 
separately, and with the total downwash being the sum of the respective downwashes. 

It is recognized that representing the effect of the rotor wake on the horizontal stabilizer 
in terms of a net upwosh or downwash is a gross simplification, but this approach makes 
the problem tractable for application to the simulation. Representing the actual flow 
field in the vicinity of the empennage with a rigorous mathematical model is not 
presently feasible. Additional dcta on pressure distributions, local flow velocities, 
etc., would be required for the development of a more rigorous model. 

Sideslip was found to reduce the magnitude of the upwash at the horizontal 
stabilizer. The downwash velocity was determined by assuming the total dynamic 
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pressure ratio to be equal to its value at zero «;ideslip. The reduction in upwash 
velocity causes a nose-up pitching moment when the aircraft is sideslipped. 

The effect of rotor wake-airframe-ground aerodynamic interference on XV-15 
handling qualities was also evaluated by observing model flying characteristics during 
semi-freeflight tests and from pilot comments during real time, pilot-in-the-loop, flight 
simulation. 

In the wind tunnel semi-freeflight tests, the model could be trimmed at the 
lowest speed achievable in the tunnel (approximate!/ 16 knots full scale). In the 16 to 
35 knots range, control was difficult but manageable. Above 35 knots, the model became 
/erv easy to trim and control, and above 50 knots the model would maintain a trimmed 
condition "hands off." The model was also flown in sideward and rearward flight at 
speeds up to 35 knots. In both cases, roll control was difficult but manageable while 
pitch and yaw control were fairly easy (even in rearward flight). The model was 
noticeably more stable in descents, including autorotation, than in level flight and 
climbs. A brief investigation of controllability at nacelle incidence angles o^ 60° 
and 30° showed the model to be easy to control about all axes. 

Real time, pilot-in-the-loop simulation was conducted on the NASA -Ames 
Flight Simulator for Advanced Aircraft. The mathematical model included the static 
instability in roll in hover, the effect of Induced Ground Effect (I GE)operotion on rotor 
performance, the wing download variation with airspeed, and the effect of the rotor's 
wakes on the horizontal stabilizer and fin. 

Although the effects of the interference were not evaluated directly, their 
influence was noted during the evaluation. Observed effects included: 

9 Increased workload to control roll when operating IGE. This effect 
of the static instability in roll was most apparent when taking off or 
landing . 
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• Stick position reversal, directional instability, and pitch coupling 
with sideslip during transition from hover to forward flight. 

e Nose down pitching when increasing collective pitch at speeds 
between 20 and 60 knots. 

All of these characteristics were similar to those observed in the scale model semi- « 

freef light tests. 

Reference 34 describes aerodynamic interference effects during low-speed 
flight close to the ground. Both flying qualities data and flow visualization results are 
presented and their implication for the design of tilt rotor aircraft is described. 


i 
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SECTION 4 


APPROACH TO THE PROBLEM 

Determination of the interference effects of the proprotor wake on the aircraft 
stability characteristics requires an appropriate model of the wake in order to analyze 
the changes in aerodynamic force included on the control surfaces by the wake. Model- 
ing of proprotor wakes is based on the vortex theory of propellers developed by Goldstein 
(Reference 35). The vortex theory is a difficult boundary value problem: given a heli- 
coidal vortex sheet, find the vorticity distribution which causes the sheet to translate 
along its axis without distortion. The shape of the vortex wake is not known in advance. 
Moreover, the stretching of the pitch of the helix has as much influence on the induced 
flow as dees the strength of the vorticity. The determination of the distortion of a vortex 
system, given the strength at the trailing edge of the blade, is a formidable task (References 
36, 37, and 38). Moreover, the vortex strength depends upon the blade loading, which 
is also unknown. This problem was treated in Reference 36; however, the computer 
program developed there was too slow for practical purposes. The approach taken in 
Reference 38 was to proceed from an assumed vortex wake geometry, including contrac- 
tion and acceleration effects. 

4.1 PROGRAM SUSSA ACTS 

The approach of Reference 38 is used in conjunction with the aerodynamic 
potential analysis technique developed by Morino (see Subsection 3,1). The technique 
is implemented in a computer program SUSSA ACTS (Steady and Unsteady Subsonic and 
Supersonic Aerodynamics for Aerospace Complex Transportation Systems), discussed in 
more detail in Appendix A. In this code the motion of the surface with respect to 
the frame of reference is assumed to consist of small oscillations with an arbitrary 
variation in time. The surface of the aircraft and its wake are divided into 
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small quadrilateral elements which ore approximated by a hyperboloidal surface defined 
by four corner points. In this process, the continuity of the surface is maintained, 
although discontinuities in the slopes are introduced. The unknown is assumed to be 
constant within each element and, therefore, the integral equation reduces to a system 
of linear algebraic equations. The problems best suited to the program are steady and 
oscillatory subsonic and supersonic (linearized, attached, potential) flows around 
simple (such as simple wings) and complex configurations (such as combinations of wing, 
fuselage, tail, and control surfaces). It should be noted, however, that the general 
formulation includes fully unsteady flow. Therefore, the problem of interference of 
the aircraft (wing-fuseiage-tail) with the rotor is a particular case of the general 
formulation presented in Reference 26. 

4.2 MODIFICATIONS TO PROGRAM SUSSA ACTS 

In order to be applicable to the aircraft-rotor interference. Program SUSSA 
ACTS requires modification in order to allow for the motion of the surface of the rotor, 
the tilt rotor aircraft geometry, and the wake geometry of the rotor. The approximate 
nature of the study on interference allows simplifying assumptions to be used for the 
aerodynamic analysis. It Is assumed, therefore, that the number of blades in the rotor 
approaches infinity, yielding an actuator disk (see Subsection 6.1). For that case, the 
aerodynamic phenomena become steady, which simplifies the aerodynamic calculations. 
The modifications made to SUSSA ACTS are discussed briefly in the following sections. 

4.2.1 ACTUATOR DISK 

The present Program SUSSA ACTS does not automate the geometry of the 
tilting proprotor aircraft rotor assembly. Hence, a preprocessor, based on data provided, 
must be developed. Because of the nature of the present study, only the rotor will be 
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considered, and this is replaced by an actuator disk. The actuator disk is approximated 
with a thin ellipsoidal configuration. A standard finite element breakup is used to 
divide this surface into small quadrilateral aerodynamic panels. The preprocessor de- 
fines the nodal numbering and its location in the reference cartesian coordinate 
system . 

On the actuator disk, the boundary condition (downwash) must be generated. 
The downwash Y on the actuator disk is related to the rotor blade downwash by 

Y = ("Uqq i + n X r) • ng/cos cr + ^ Cj QR tan o (1) 

where 

a = argle of attack of the blade 
ng = normal to the blade 

n = angular velocity of the rotor 
Uqq = free stream velocity 

Cj = rotor thrust coefficient 

R = rotor radius 

In obtaining this inflow, the components of tho perturbation velocity tangent to the 
disk are assumed small. Details of this formulation are presented in Subsection 6.2. 


4.2.2 TILT ROTOR AIRCRAFT PREPROCESSOR 

The original preprocessor of SUSSA ACTS was designed for conventional 
wing-body-tail configurations. The wing was assumed to be located at mid fuselage, 
and the profile of the fuselage was assumed to be circular. In addition, the empennage 
was limited to only one vertical tail. Consequently, the geometry of the nose, fuse- 
lage, and empennage was modified appropriately to model the tilt rotor aircraft. 

These modifications affect the final nodal breakup of the tilt rotor aircraft. 
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4.2.3 ROTOR/AIRCRAF WAKE GEOMETRY 

The surface of the wake emanating from the trailing edge of the wing and tail 
Is assumed to be known. These wakes are treated as straight vortex lines from the 
trailing edges of the wings and tail and parallel to the direction of the flow. The Kutta 1 

condition is imposed at the trailing edges. 

The wake of the actuator disk (rotor) is composed of infinite vortices emanating 
from the rim of the disk and one additional vortex emanating from the center of the disk. 

The control vortex represents the vorticity shed by the blade near the blade root. Details 
of the wake geometry of the actuator disk are given in Section 6. 

• Two-Vortex Model (Helicopter and Initial Conversion Modes) 

When the rotor has a small incidence angle with respect to the free stream 
(less than 30°), the wake rapidly rolls up into two concentrated vortices similar to the 
tip vortices of a monoplane wing (References 39 and 40). 

On this basis, the wake model chosen for the helicopter mode of tilt rotor 

I 

flight consists of two vortices trailed from the ends of the lateral diameter of the rotor 
disk, as shown in Figure 6(a). The strength of the vortices is obtained by assuming 
that the rotor behaves as an elliptical monoplane wing of span 2R (the Glauert hypothesis) 
whose lift, i.e., thrust, T, is given by 



where 

p = air density 
V = flight velocity 

= mid-span value of elliptical circulation distribution 

y = spanwise coordinate measured from mid-span 
R = rotor radius 
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(b) Resultant Flow Velocity 


Figure 6. Two-Vortex Model. 
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Making the transformation y = R cos 0 , dy = -R sin 9 d 9 , the thrust becomes 


T 


r,„ 

Jn 


pVr^R I sm 0d0 


n 


2 


pvr R 

o 


( 3 ) 


Then the strength of the trailing vortices is given by the strength of t^e mid-span 
circulation, i.e.. 


- pVR 

2 


In the far wake the two vortices are assumed to be convected downward at 
one-half the local downwash velocity, due to vortex sheet roll-up effects, as suggested 
by the results on page 45 of Reference 39. The resultant flow velocity is then given by 
the vector (see Figure 6(b)) 

V' = [(V + 2v cosGfj^), (vsinoj^)] (5) 


where, from Reference 41, page 185, the induced velocity, v, is 

c^nR 


and 


V sin Qf|^ 

M = 

HR 


V cos cy|^ + V 

nR 


( 6 ) 

(7) 

( 8 ) 
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I 2 2 

pnR-'inRr 

This model is assumed to be valid for the following conditions: 

V = 0 to 130 knots 

= 60 to 1 00 degrees 
K 

• Vortex Cylinder Model (Final Conversion and Cruise Modes) 

When the rotor has a large incidence angle with respect to the free stream 
(more than C0°), the wake tends to assume the shape of a skewed cylinder of vor- 
ticity (Reference 39). 

On this basis, the wake model chosen for the final conversion mode of tilt 
rotor flight consists of a skewed cylindrical surface of vorticity trailed from the edge 
of the rotor disk, as shown in Figure 7. The steady-state vorticity distribution is 
obtained from the expression derived in Reference 42 for the bound vortex strength of 
one rotor blade: 


r 

0 


2T 

pOR^b 


(10) 


where 

Q = rotor rotational speed 
b = number of rotor blades 


Then the vorticity distribution on the surface of the cylinder is given by 


£ 

dx 


br 

o 

Ax 


( 11 ) 


where Ax = length of vortex cylinder generated during one revolution of the rotor, 
and measured along the axis of the cylinder. 
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Also, it is noted that 


Lx = V" Lt 


( 12 ) 


where V" denotes the mean convection velocity of the far wake and At = for one 
rotor revolution. 


Therefore, Equation (1 1 ) becomes 


IT 

dx 


br n 

o 

2nV" 

T 

priR^V" 


(13) 


The convection velocity of the portion of the vortex cylinder in the far wake 
is given by the vector (see Figure 7). 


V" - [ (V + 2v cos o-j^), (2v sin ofj^) ] 


(14) 


where v is given by Equation (6). 

This model is assumed to be valid for the following conditions: 
V = 13C Units to V 

max 

ccj^ - 0 to 60 degrees 
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SECTION 5 


INCOMPRESSIBLE POTENTIAL AERODYNAMICS FOR 
WING-FUSELAGE-TAIL-ROTOR WAKE INTERFERENCE 

In this section, the general formulation for incompressible potential aerody- 
namics for separoted flows is presented. It may be worth noting that for incompressible 
fluids potential flows are obtained if the fluid is inviscid and the velocity flow field is 
irrotational at time t = 0. Accordingly, these hypotheses are invoked here. Since the 
analysis of rotors was not included in previous descriptions of the formulation for SUSSA 
ACTS (References 43 through 47), a new formulation is presented here, starting from 
basic principles. 

5.1 INCOMPRESSIBLE INVISCID FLOWS 

★ 

Assuming the fluid to be [1] inviscid, and [2] incompressible, the motion is 
governed by the Euler equations 


— + V- V^V = -lVp 


and the continuity equation (for incompressible fluid) 


(15) 


V • V = 0 (16) 

where V is the velocity vector with respect to a prescribed frame of reference, p is 
the pressure, p the (constant) density, t the time, whereas 


V - ~ ^ ,“^,7 ^ 

V = I + I + k — 


ax 


ay 


a z 


(17) 


is the del operator. Equations (15) and (16) form a system of four partial differential 

equations for tour unknowns V , V , V , and p . 

X y z 


*Numbers in brackets indicate new hypotheses introduced in the formulation. 
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In order to complete the Formulation of the problem, the boundary conditions 
at infinity on the body and on the wake must be obt< vl. Here it will be assumed that 
[3] the unperturbed flow consists of uniform translation, say in the x-axis direction. 
Hence, the boundary condition at infinity may be written as 

9 = U^T (18) 

In particular = 0 for helicopters in hover. 

On the body it is assumed that (4) the surface of the body is impermeable, 
i.e., the norma! components of the velocity, V, of the fluid, and of the velocity, Vg, 
of the body coincide on a point, P, on the surface. Eg, of the body: 

(V - Vg) • n = 0 (for P on Eg) (19) 

where n is the normal to Eg at P . 

Next consider the boundary condition on the wake for lifting bodies. Here 
it is assumed that [5] the body has certain separation lines (in particular sharp trailing 
edgesl and that the velocity flow field has surfaces of discontinuity (wakes) emanating 
from these separation lines (or trailing edges). Indicating with 1 and 2 the two sides 
of the wake, let n be the outward normal on side 1 and let 

A f = f, - f2 (20) 

indicate the discontinuity of any functions, f, across the wake surface. (For the classi- 
cal wing-wake, 1 and 2 correspond to upper and lower sides, respectively; n is the upper 

normal and Af = f - fi ). The boundary condition on the wake-surface, E..,, 
upper lower ' w 

is that there is no pressure discontinuity on a point, P, on the wake 

Ap = 0 (for P on E^) (21) 

This condition can be expressed mare conveniently in terms of the velocity 
discontinuity; however, this requires the introduction of the concept of Lamb surfaces. 

-30 - 


ABROaPACB ■VSTBMa, INC. • one vinc BmoK park • auRLiNOTON, maseachubbtts □teas • (bi7) a7a- 







These surfaces may be avoided in the case of potential flow and, therefore, the 
boundary condition is given here for potential flows only. 

Equations (15), (16), (18), (19), and (21) completely define the problem in 

terms of four partial differential equations for four unknowns V , V , V and p, with 

X y z 

the corresponding boundary conditions. 

5.2 DIFFERENTIAL FORMULATION FOR INCOMPRESSIBLE POTENTIAL FLOWS 

In the preceding section, the flow was assumed to be [1 ] incotr^ressible, and 
[2] inviscid. Next, assume [6] the flow to be initially irrotationol: then the flow re- 
mains irrotational at all times and, therefore, there exists a potential function, $, such that 

V = V$ (22) 

In this case, Equation (15) may be integrated to yield Bernoulli's theorem* 


+ L |7*p + B. = 1 ij2 + 
at 2 ' I p 2 ® 


00 

p 


(23) 


Tlie constant on the right hand side is obtained from the conditions at infini:/. 
Equation (18). 

Furthermore, Equation (16) may be rewritten as 


= 0 


(24) 


*The x-derivative of Equation (23) yields 


i |v*|2 + £ 

3x \ at 2 "I I p 


^ a^ a $ ^ 1 ap 

az axaz p ax 


at 


^ a^ ^ ^ a^$ 


ax g^2 dy axay 


= ( V- vw 

at 


+ 1 ^ = 
p ax 


in agreement with the first component of Equation (15). 
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f 




2 

where V Is the Loplacian operator. 

Similarly, Equation (18) is equivalent to 


# = X (for P at oo) 


while Equation (19) may be rewritten as 


3# v> 

_ = • n 

3n ° 


(25) 


(26) 


Next consider the boundary condition on the wake. Equation (21). Combin- 
ing Equations (21) and (23) yields 

1 - ~ ^ 2 ) ~ ^^*1 ‘^^1 ■ ^*2 ' ^* 2 ) " ° 

3t '2 


or 




A*= 0 


(28) 


where 


A# “ " ^2 


(29) 


whereas 


,A= 1 (v*, L (7, .Vj) 


(30) 


is the average value of the velocity: this average value may be attributed to the point, 
P^, on the wake. Equotion (28) may thus be inlarpreted as 

A$ (P ) = constant in time (31) 


w' 


where P^ is the physical point of the wake being convected with velocity V^. Note that 
Equation (28) implies that the wake is tangent to In particular, in the steady state. 
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the wake is composed of the streamlines of the points P emanating from the trailing edges 
or other separation lines. Equation (24) is a partial differential equation for i with bouna- 
ary conditions given by Equations (25), (26), and (28). Once i is known the velocity is 
given by Equation (22) while the pressure is obtained from Equation (23). It may be worth 
noting that the potential aerodynamics formulation considered in this section is consider- 
ably simpler than the formulation based upon physical variables (V and p) since the former 
gives one linear partial differential equation instead of four nonlinear ones. 

5.3 PERTURBATION POTENTIAL 

In order to use the Green theorem approach, it is convenient to have homo- 
geneous boundary conditions at infinity. This can be accomplished by introducing the 
perturbation potential 9 , related to the potential $ by 

* = X + cp (32) 

Then, Equotion (24) yields 

= 0 (33) 

while the condition at infinity. Equation (25), becomes homogeneous 

cp = 0 (for P at 00 ) (34) 

On the other hand, combining Equations (26) and (32), the boundary condi- 
tion on the body is now given by 

= (-Uj + ■ n (35) 

on 


whereas, on the wake (see Equations (28) and (31)) 




J Acp = 
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or 

A cp (P^) = constant in time (37) 

Furthermore, once cp is known, the perturbation velocity 

V = Vcp (38) * 

may be evaluated and then the velocity is given by 

V = T + V (39) 

Accordingly, the Bernoulli theorem may be rewritten as 

^ 4 u ^ 4 1 |Vcp|2 + P = ^ (40) 

9 t °° 9x 2 p p 

It may be worth noting that the above equations could be obtained in a 
slightly different way, i.e., by assuming the frame of reference connected with the 
undisturbed air: the velocity of the body in such a frame of reference would be given « 

by -U i + Vn, in agreement with Equation (35). 

00 D ^ 

5.4 GREEN'S FORMULA FOR LAPLACE'S EQUATION 

Consider two orbitrary functions f and g and note that 

V • (fvg) = Vf • Vg + fv^g 

_ _ _ _ 2 

V • (gv 0 = vg • v f + ov f (4i) 

or 

V • (fvg - gvO = f v^g - gv^f (42) 

According to the Gauss th'sorem, for any arbitrary vector a, t 


-34 - 


ABRO'BPACB SVaTBMai INC. • one vine brook park • burlinqton, Massachusetts o'leoa • (si7) a7a-7Bi7 





IJJ ^ • “■“''i = # a • n ds (43) 

V. ' I 

^ where I! surrounds V. and n is the inward directed normal (Figure 8). Combining 

Equations (42) and (43) yields 



Figure 8. Geometry for Green's Formula for Vj Internal to S . 

Ill 9 f) d V = 0(f vg - g V f) • n dr 



s 

Next assume that both f and g satisfy the Laplace equation, i.e., 

V^f = V^g = 0 (inV.) (45) 

Then Equation (44) yields 
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( 46 ) 




which is the desired Green's formula for the Laplace equation. 

Next choose the function g to be the simplest function which satisfies the 
Laplace equation: if g is only a function of the distance, r, from a specified point, 
P*, then the only such function that satisifes the Laplace equation is 


g = - + B 
r 


(47) 


where 

r = I P - P* I (48) 

Note that P* must be outside the volume, V»; otherwise Equation (45) is not satisfied 
at P* . Note also that the value of B does not affect Equation (46) under the condition 



which is obtained from Equation (46) with g = 1. Therefore, the convenient value 
B = 0 is used here. For this case the function g represents a source. It is convenient 
to choose A such that g represents a unit source; that is, a source with flux equal to 
one. This is obtained with A = “1/4 tt, i.e.,* 


''For, in this case, the flux through a spherical surface of radius R 


ffy -n dr ^ ff ^ 




r = R 


4n R 


, is given by 
4nR^ = 1 
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g = - 


4n I 


(50) 


Combining Equations (46) and (50) one obtains the desired Green formula for 
the Laplace equation 


(£)-^ 


(51) 


5.5 GREEN'S THEOREM FOR INCOMPRESSIBLE 
POTENTIAL AIRODYNAMICS 

In order to obtain the Green theorem for inconpressible aerodynamics 
(Laplace's equation for the exterior problem, i.e., for Vj outside the surface J^q)/ 
it is convenient to consider the surface E composed of three branches as indicated 
in Figure 9. Note that the inward normal for V. Is outward to Eq . E| is o 
spherical surface of rad' ■ .ind center P* , while E 2 is a spherical surface of 
radius R_ and center P* . 

Then Equation (51), for f = cp, reduces to 


^ |_ Bn \ 4nr / Bn \4rrr/ 


Eq + E] + E2 


Note that as the radius e of E^ goes to zero 
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Figure 9. Geometry for Green's Theorem for Potential 
Aerodynamics of Incompressible Flows. 
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i- 



f, 

t' 

\ 

?' 

c 

t 

1^ 

? 


! im 

G-»0 E 


dE 


= lim 
e-* 0 

= lim 
e~* 0 

= lim 
e-*0 


/^\ -!£ (^) 

^ 1_ 9n \4TTr/ dn \4TTr/ 


’•#- (—) 

^ Sr \ 4tt r / \Sn/* ^ 4rrr 


L E 


4tt e 


(ftds -i- (ft dr' 

^ \ Sn /* 4rt e ^ 

^1 , 


9* + (— ) e ^9 


whereas, as the radius goes to infinity 

^2 


'1 


(53) 


lim 
-* cx> 


So ^2 


(ftLA (^).^ 

^ j_ Sn \ 4it r / Sn \ 4n r / _ 


dE 


lim 




9 ■*■ 

^00 




R* 


= 0 

under the conditions^ 


(54) 


lim cp - 0 

R^ -»Q0 
^2 


(55) 


f 


^It is verified a posteriori (Equations (67) and (68)) that these conditions are 
satisfied. 
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. M , 



r/»Sf 


lim ^ Rj, = 0 


iL -»CD 
^2 


9n 


(56) 


Therefore, as e tends to zero ond R_ tends to infinity. Equation (52), 

^2 

according to Equations (53) and (54), yields 


- ■«[!:(£)-•£ (a] 


dSf 


(57) 


It may be worth noting that if the point P* is inside £q, then there is no 
need for the surface and, therefore, in this case Equation (52) yields simply 


0 = 


# [g (ir;) ■ • ^ bz) ] "'o 


(58) 


Equations (57) and (58) may be combined by writing 


*■' ■#[;• (,T.) -V. (.T,) ]"• 


(59) 


E* = E(P*) = 




(60) 


(P* outside Eij) 

( 0 (P* inside E^) 

Next let Eq be identified with the surface surrounding the fuselage, rotor 
and wakes emanating from the trailing edge of wing, tail and rotor blades. Then 
Equation (59) may be rewritten as^ 


4-.E..P. = - (ft ris i - /i\ 

v'It' r 9n \ r / 




dE 


1/ .. i (1) 


dE 


(61) 




^Note that the wake surfaces Ey^, e| and Eo are open surfaces. The source integrals 
on the wake surfaces are identically equal to zero since the values of Scp/dn on the 
opposite sides of these surfaces are equal to zero. 
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where ZIg is the (closed) surface of the body (wing-fuselage-tail), is the (closed) 
surface of the rotor blades, is the (open) surface of the wing wake, Zj is the (open) 
surface of the tail wake and Z^ represents the wakes of the rotor blades. Furthermore 
(see Equation (20)), 


A cp = cp^ - cp2 


(62) 


while n is the normal on the side 1 of the wake. Note that Acp is known from Equation (37). 

In addit’on if P.. is on Z„ + Z„ then a. shown in Appendix C of Reference 45, 

D K 

E* - 1/2* Thus Equation (60) may be rewritten as 

1 


E* = E(P*) = { 1/2 
0 


(P* outside Zg + Zj^) 

(P* ^R) 

(P, inside Zg+Z|^) (63) 

Note that on Zg and Zj^, Scp/Sn is prescribed by the boundary conditions, Equation (35). 
Thus, for P* on Z , Equation (59) is an integral equation relating cp to the prescribed 
hi /Sn . This equation is used to analyze the potential flow for interaction of fuselage 
and rotor. 


5.6 CONDITIONS AT INFINITY 

In this section Equations (55) and (56) are verified. Note that R^ goes 
to infinity, .Equation (57) yields 


lim cp* = lim 


"r."“ E 


^ |_3n \^trry 9n (4rrr 


I 


dZ 


lim 


R^ -00 
^2 


4nr ^ 3n 
- ^0 


dZ - 


1 




dZ 


41 - 


4rrr Jq 
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since the flux through Eq is equal to zero, i.e.. 



while 



dZ 


finite 


Equation (61) implies 


( 65 ) 


( 66 ) 


cp = — + 0(r ^) at CO 

r 


and 


(67) 


^ - il + 0(r"^) at 00 (68) 

5n ^3 

in agreement with Equation (56). 
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SECTION 6 

INCOMPRtSSIBLE POTENTIAL AERODYNAMICS FOR COMPLETE 
CONFIGURATION WITH ACTUATOR DISK 

The formulation presented in the preceding section is general and may be ap- 
plied, in particular, to rotors, it may be noted, however, that the rotor-fuselage inter- 
action is an unsteady periodic phenomenon. Periodicity is obtained when the rotor moves 
through an angle equal to 2n divided by the number of blades, b. If b goes to infinity, 
the phenomenon becomes steady. Therefore, given the preliminary nature of the present 
analysis, it is convenient to consider the limit of the rotor fuselage analysis as the num- 
ber of blades goes to infinity. The model representing an infinite-blade rotor is called 
an actuator disk. The modifications necessary to extend the formulation to include actu- 
ator disks are presented in this section. For simplicity an actuator disk by itself is con- 
sidered first. The extension to the complete configuration is considered in Subsection 6.3. 

6.1 ACTUATOR DISK 


Consider a rotor composed of b blades. Let b go to infinity, while the chord 
of the blades goes to zero. As the aspect ratio of each blade goes to infinity, the 
voiiicity shed by each blade is more and more concentrated near the root and the tip of 
the blade. In the limit, only tip vortices and one central vortex exists. Therefore, the 
actuator disk may be represented as a series of infinite biades with a wake composed of 
infinite vortices emanating from the tip of the blade and one additional vortex emanat- 
ing from the center of the disk. Thus, the flow field is potential outside the surface Zq 
shown in Figure 10. Therefore, Equation (59) is applicable in this case, and yields 




^OOli quality 


"" C i (i) " 


(69) 


D B 
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Figure 10. Surface Eq for Actuator Disk Analysis. 

The Surface Surrounds the Actuator Disk, 
the Central Vortex and the Shed Tip 
Vortices. 




where Ep is a (closed) surface surrounding the blades of the rotor and representing the 
actuator disk, while is the surface of the disk wake emanating from the circle 
described by the rotor-blade tips and Eg is a branch surface connecting the central 
vortex to Ep . It may be worth noting that Equation (69) is formally equal to 
Equation (61). Therefore, the method presented in the preceding section is applicable 
to actuator disks as well. The two differences are that 9cp/9n cannot be obtained from 
the usual boundary conditions on the surface of the body and that A cp on Eg cannot be 
determined. This last point is examined in Appendix B. The new boundary conditions 
for actuator disk aerodynamics are presented in the next subsection. 


6.2 ACTUATOR DISK BOUNDARY CONDITIONS 

The boundary condition on the rotor blade according to Equation (35) with 

Vg = n X 7 (70) 

is 

it , (-u<„i + nx;) . nj (71) 

3ng 

where ng is the normol to the blade. On the other hand, the boundary conditions used 
on the surface of the actuator disk ore expressed in terms of the normal wash, Sip/Sn^ , 
on the surface of the disk (see Figure 11). Therefore, the relationship between Bcp/Sng 
and 3cp/^npis necessary and is obtained here. Let Sg and s^^ be the directions of the 
components of the velocity tangential to the surface of the blade and the disk, respec- 
tively. Then 


V cp 




7 ^ ^ 

'■^B ^*B 




(72) 


- 45 - ORIGINAL PAGE : ' 

OF POOR QUALm 

ABROSPACE BYBTBMB, INC. • onb vine brook park • buruinoton, MASSACHuaeTTB oigo 3 • (317; a?' 




Figure 11. Boundary Conditions for the Actuator Disk. Eq is 
the Surface Surrounding the Actuator Disk. 






Note that 


^2 vcp • UrN - n„ • s„ • n 


9n, 


'D 


9n, 




(73) 


and 


^ ng • Sd + ^ sg . Sp 

3sd 3ng 5sg 


(74) 


Eliminating 9 cp/bsg between Equations (73) and (74) yields 

9cp ^ ^ ^'D • ~ ‘ ^ 


9np ang 


* *D 




nB X Sb . np X s^ ^ ^ * ''d 9cp 


®B * ®D 


9nc 


*B ' *D 


(75) 


Next note that 


iSL / = tan P 

9 sp/ 3 np 


(76) 


where 6 is the angle between 79 and np. Therefore Equation (75) may be rewritten as 


"D \ ®B ’ *D / ®B ’ ®D 


3 


3ng 


(77- 


which is the desired relationship between (necessory for the solution of the problem) 

3 

and — 2- (known from the boundary conditions on the blade. Equation (71). 

3ng 
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Next consider a simpler expression for Equation (77): note that usually the 
radial component of is negligible. Then the configuration on a cylindrical surface 
is represented in Figure 11, In this case 


ngxsg 



= cos O' 


= sin a 


(78) 


and Equation (77) reduces to 

Sep _ 1 1 3cp 

Qn^ cos a 1 - tan Of tan p Sn^ 


(79) 


where a is the angle of the blade section for a given angular position.* 

The boundary condition on the actuator disk can be written in terms of rotor 
parameters. It is shown in Figure 11 that 

V 


tan P = 


Sep 

Sn« 


(80) 


with denoting the circumferential velocity component. The closed line integral 
of the velocity vector for the actuator disk yields an additional relationship (assuming 
N to be independent of the circumferential location). 


Aepg = ^ 


• d s == ?rr R V 


(81) 


*|t is emphasized that different values of f in the boundary condition do not affect the 
lift (thrust) on the rotor, although the flow through the disk is changed. The lift 
(thrust) on the rotor is dependent upon the value of i.e., the discontinuity of 

the branch wake. 
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where represents the circumferential discontinuity on the branch wake. Combining 
the above two equations yields 


tan p 


% 

2tt R 


1 


f3nr 


'D 


Hence, Equation (79) can be reduced to 


(82) 


acp 


1 


9n, 


“Pb 

2nR 


tan a 


(83) 


•q cos or 

Under the assumptions used in Appendix C, the discontinuity, Acpg, is proportional to 
the thrust, T (see Equation C-16). 

T 


Acp 


B 


1 pQR2 

2 


(84) 


Therefore, the actuator disk boundary condition becomes 

^ + - T- . (85) 

or (see Equation C-17) 

^ = — 1 - ^ - i C-|.nR tana ( 86 ) 

9 np cos a gn^ 2 

This is the desired expression which relates the boundary condition on the actuator 
disk to the rotor thrust and the downwash on the rotor blades. 
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Using classical blade element theory (Reference 41), the thrust coefficient 
(for rotor with large inflow) can be obtained in terms of specific rotor parameters. 

The exact relationship is given in Eq- Jtion (108). In obtaining this expression, the 
profile drag is neglected, the angle of attack of the blade clement is represented by 
its sine, and the pitch of a wisted blade is approximated by its value at three-quarter 
radius. 

6.3 FORMULATION FOR COMPLETE CONFIGURATION 

Consider a surface, 2 , surrounding the aircraft and its wakes as well as the 
actuator disk and its wake. Outside the surface 2 the flow is potential and, therefore. 
Equation (61) may be applied, with the surfaces of the disk , and of its wake, 

Zp replacing the surfaces of the rotor, and its wake, Fj^. Therefore, Equation (61) 
may be rewritten as 


4ttE*cp* - (ft ^ 1 - cp — 

^ an - 9n V > 




dZ 


If “"vO 

y' + y- °n r 

T D 


dZ 


(87) 


The value of E* is given by the equivalent of Equation (63), i.e.. 


1 for P* outside Z_ + Z-. 

B D 

E* = E(PJ = <; 1/2 for P* on Zg + Zjj 

0 for P* inside Zg + Z|^ 


( 88 ) 


- 50 - 


ABROBPACe BVSTBMB, INC. • ONE VINE BROOK PARK • BURLINGTON, MA8BACHUBETTB 01S03 • (ai7) 



Note that on Zg, Scp/Sn is prescribed by Equation (35), whereas on 
Scp/Bn = c'cp/dnQ is given by Equation (79) with 9cp/Bng given by Equation (71). 
Note also, that with the introduction of the steady state approximation accomplished 
through the use of the actuator disk. Equation (37) simplifies considerably to yield 

A 9 (P ) = constant along streamlines (89) 

w 

In addition, the geometry of the wake may be estimated as illustrated in 
Subsection 4,2.3. Therefore, on Zg + Zp, Equation (87) is an integral equation 
relating the unknown values of cp on the surfaces Zg and Zp to the corresponding 
value of the known normal derivative Bcp/9n (normal wash). The numerical solution 
of the integral equation is considered in Sev':tion 7. 
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SECTION 7 

NUMERICAL FORMULATION 


The ground formulation for incompressible potential aerodynamics was pre- 
sented in Section 5. While the formulation Is valid for fully unsteady aerodynamics, 
considerable computer time saving is obtained by the use of steady-state aerodynamics, 
which is achieved by replacing the rotor with an actuator disk (Subsection 6.1). The 
integral equation for the complete configuration is given in Subsection 6.3. The numeri- 
cal formulation for the solution of the integral equation is outlined in this section. 

The surfaces of the aircraft, the actuator disk, and their wakes are divided 
into small quadrilateral elements. Each element is replaced by a hyperboloidal surface 
defined by the four corner points of the element. In this process, the continuity of the 
surface is maintained. The unknown is assumed to be constant within each element and, 
therefore, the Integral equation is approximated by a system of algebraic equations. 


7.1 APPROXIMATE ALGEBRAIC SYSTEM 

The results of Section 6 may be summarized by saying that the problem of 
incompressible potential aerodynamics for the complete actuator disk -aircraft configura- 
tion is governed by the integral equation obtained by combining Equations (87) and (88) 
on the body surface to obtain 


2ncp* = - 






+ 


II 


e' +s' +r' +s' 
T D B 



dS 


(90) 


fAGt KAMi 
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In order fo obtain an opproximote solution for the integral equation, the 
surfaces Tg, Eq, Ey, and Ep are divided into small quadrilateral elements, E|^ . 
The values cf tp and dcp/9n are assumed to be constant within each element. The 
collocation method is then used, that is. Equation (90) is satisfied at the centroids, 

Pj^, of the element, E^ . This yields 



‘^hk "'''hk^ I'Pki “ ^^hk^ 



+ Acp^W® 


(91) 


where 5|^|_ is the Kronecker delta, 



(92) 


(93) 


In addition, w^j^ = 0 for the elements not in contact with the trailing edge, while for 
the elements in contact with the trailing edge. 


hk " 



(94) 


where Ej^ is the strip of the wake (bounded by two streamlines) emanating from the 
element Ej^. The upper (lower) sign must be used for the upper (lower) side of the 
wake. (In deriving Equation (91), the value of cpyy is approximated with the value 
of cp at the centroid of the elements adjacent to the trailing edge. This is reasonable 
in view of the Kutta condition.) Finally, AcpgWj^ with 
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= 


2n 


ff - 

5n 




dS' 


(95) 


■JP. 


is the contribution of the branch wake. Note that Aep = Zicpg is constant on the 
branch wake and, therefore, there is no need to divide Eg into small elements. Note 
that Acpg cannot be determined (see Appendix B), but has to be prescribed. Therefore, 
it is left on the right-hand side of Equation (91). 

As shown in Appendix B, Equation (91) has the determinant equal to zero; 
to be more precise, one equation is a linear construction of ^he other. Therefore, in 
order to find the solution, one equation is dropped and one unknown is set equel to 
zero. Then, Equetion (91) can be solved nuinericolly to yield the remaining values of 
the unknowns, cp|^ . In all the results presented here, the pressure is obtained from the 
Bernoulli theorem, Equaticn (40). 


7.2 


HYPERBOLOIDAL ELEMENT 


The surface elements E|^ are approximated by a hyperbolic paraboloid 
(hyperboloidal element. Figure 12) 

P = Pq + ?Pi + ^P2 ^ - ^P3 


(96) 


where p. are obtained In terms of the locations of the four corner points as 


Po' 


"1111 

1 ' 

Pi 

_ 1 

1 1-1-1 

|P+- 

^ - 

^ " 4 

1-1 1-1 


1 P2 ' 

1 

1 1 

- 


1 -1-1 1 

1 _ 

V P 3 

' 

\ P__ / 


(97) 
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Nole that quadrilateral hyperboloidal elements can be combined to yield a closed sur- 
foce. Using these elements, the coefficients Cj^j^ ond bj^j^ con hi evaluated analytically. 


as 




D' 


( 98 ) 


‘"hk = 's(i/i)-'s(^' -^)-isH, I' + lsH, -1) 


with [using -tt/2 < tan~^ ( ) ^ ^/2] 


lQ(^/‘n) = (l/2n) ton V(rxa^ •rxa 2 )/(| r [ 7 • a^ xa 2 )j 


( 99 ) 


Ooo) 


and (neglecting Sn/9| and 5n/oTl) 


:(?,T1) 


2tt ( 


-r X a, • n sinh 


l'I±] 


rxa 


. - - 1 • L-l| ■ °2 

+ rxa^ •*n smh I 

z I I 


I r X a^ 


where 


+ 


r> tan 

P 


X 

1 

°1 

• 7 x ^2 ) 

1 

\|r 

Ir 

• °i X 02 y 

1 


'(?, '0) = ^ " Pq ^P] ^P2 *^P3 "^h 

O] (?,'fl) = 3r/S? = + Tlp^ 

02 (5'^) = 9 it's T) = ?2 + 5P3 
n (|,T1) = (a^ X 02 )/| a^ x C 2 | 


similar expressions are used for Wj_|^ 
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SECTION 8 


DETERMINATION OF FLIGHT CONDITIONS 

The effects of the aerodynamic interference of the rotor wake and the air- 
craft on the stability and control characteristics of the tilt rotor aircraft art expected 
to be important. In par»-icular, the effect on the longitudinal trim of the aircraft, 
especially through the transition and conversion modes, may be significant. In .he 
trim condition, the total longitudirKil farces and moments acting on the aircraft must 
be in equilibrium, that is, they must sum to zero. For a given mast angle, flap setting, 
and stabilizer incidence, the aircraft is trimmed by appropriate displacement of the 
longitudinal and collective pitch controls. 

Because of the complexity of the forces and rr’oments on proprotor aircraft, 
an Iterative solution is usually required. In which the displacements of the controls 
required for trim are first estimated, and then Incrementally changed to achieve force 
and moment equilibrium of the aircraft within some acceptable range, soy ± 50 lb or 
ft-lb from zero. However, for the purpose of the present study, an alternate method 
of determining the trim conditions in the presence of aerodynamic interference was 
selected. This method was based on incremental chonges from a standard case in 
order to estimate the values of collective pitch and rotor incidence for other flight 
coses of particular interest. The standard case was based on wind tunnel data of a tilt 
rotor aircraft in the trimmed condition. 

Estimates of collective pitch and rotor incidence under various conditions 
of weight, altitude and flight path angle are made through an incremental anolysls 
based on this standard case (Subsection 8.2). This involves calculating the increments 
in vehicle drag and rotor thrust coe.^icienJ based on an assumed vehicle drag polar, and 
the rotor blade modr ' . The rotor blade model is based on the application of classical 

«acsDiNa fAGt ujm mnnm» 
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blade elenTient theory to compute the thrust and flopping motion of o rotor with large 
inflow (Subsection 8.1). 

8.1 APPROXIMATE ANALYflS OF THE THRUST AND FLAPPING 
MOTION OF ROTORS '.v'lTH LARGE INFLOW 

Applying classical blade element theory to the section of rotor blade shown 
in Figure 13, neglecting the effect of profile drag, and assuming that the angle of 
attack of the blode element con be represented by its sine, and that the pitch of o 


ZERO LIFT LINE 



Figure 13. Typical Blade Element at Radius, 0.75 R. 


twisted blade con be represented by its value at three-quarter radius, the rotor thrust 
perpendicular to the tip path plane can be expressed as: 


where 


T 


_b 

2tt 



P 1.2 


- U acsin [ 0 q 75 “0l cos 0 drdt 


(106) 


b 

P 

U 


u_ 


u. 


number of blades 
air density 

resultant air velocity 
\nR + pORP^ cos (f 
Or -t- p QR sin t 


u2 . 
L p 



1/2 
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\ = 


U sin i + V 
inflow rotio = — 


= advance rafio = 


nR 

U cos i 


00 


nR 


I 

V 


rotor incidence angle 

rotor induced velocity (assumed constant) 


Cy = rotor thrust coefficient = 


pttR^ (nR)^ 


2 [,2 . , 2 - 


1/2 


n = 

R = 

^ = 
r = 
a = 
c = 


0, 


0.75 


rotor rotational ^eed 

rotor radius 

blade coning angle 

radius to blade element 
blade lift curve slope 

blade chord, assumed constant 
blade pitch ct three-quarter radius 


= 0 , 


0.75 " "1 


+ 9, cos t + 0^ sin i|r 


0 


0.75 

c 

s 

t 


C s 

= blade collective pitch at three-quarter radius 
= blade lateral cyclic pitch 

= blade longitudinal cyclic pitch 
= blade azimuth angle 


-1 -1 

0 = blade inflow angle = sin ' — P or cos” - 

U U 


Expanding and substituting for cos 0 and sin 0 , 

.2tt R 


b r ^ D r ? 1 

T = — I j — ac U-r sin - U Ut cos 0,« drd^ 
2n L I 2 L T 0.75 p T 0.75J 


(107) 
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2 2 . 

Nondimensionalizing the thrust expression by pnR (nR) , and integrating, the thrust 
coefficient is obtained: 


C 


T 


_ aa 

p'" ®0.75 I 

1+^ 2\ 
1 + -P 

2 

3 I 

( 2 I 




+ (cos Sq 75 + X sin 6|J 75) 


(108) 


where the solidity, a, of o rotor is defined as the ratio of the total blade area tc the 
rotor didc area. For rectangular blades. 


bcR _ ^ 

"I 


(109) 


All quantities are referred to the tip path plane. 

Applying blade element theory in die same manner as for thrust, and assuming 
negligible spring restraint about the flapping hinge, the equation for equilibrium of a 
single blade about the flapping hinge is 





( 110 ) 


where 1^ = blode moment of inertia about flapping hinge, 

2 

Nondimensionalizing the flapping expression by l^n , integrating, and equating like 
harmonics, the equilibrium blade coning angle and cyclic pitch with respect to the tip 
path plane ere obtained: 


e 


o 


Y 

7 


^ 0.75 

4 


P-01 

(1 +P^)+ ^ 

3 


cos 0 


0.75 


, 3 , . . 
+ — X sin 0 

4 



X 

3 


cos 00^75 


( 111 ) 
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V3 

1 + T 

= P [^°»Qq,75 - T 0 

’s 1 + y + y »on 0Q 75 


( 112 ) 


(113) 


in determining the equilibrium of the aircraft as a whole, It is necessary to 
determine the flapping of the blades with respect to the hub plane. It can be shown by 
an analysis similar to that described above, that the following relationships exist: 


where 


1 


'HP 


1 + 


Pi 


3 ^ ^0.75 *HP 

w 1 2 

1 + y p 


®1 

*HP 


1 


1 + 


P, =01 


3 ^ ^°”^0.75 ^HP 
1 + 2 M 


014) 


(115) 


Pi = blade lateral flapping angle with respect to hub plane 

*HP 

01 = blade lateral cyclic pitch angle with respect to hub plane 

?! = blade longitudinal flapping angle with respect to hub plane 

*^HP 

01 = blade longitudinal cyclic pitch angle v/ith respect to 

S|_|p hub plane 


and all other quantities are referred to the tip path plane. 
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8.2 


SELECTION OF FLIGHT CASES FOR ANALYSIS 
INCLUDING ESTIMATED VALUES OF COLLECTIVE 
PITCH AND ROTOR INCIDENCE FOR TRIM 


In studying aerodynamic interference effects on tilt rotor aircraft, it is first 
necessary to establish trim condi'-ions for the various flight coses of interest. Parameters 
to bo varied ore gross weight, altitude, climb angle, and center-of-gravity position. 

During flight, o typical proprotor aircraft has three distinctly different modes 
of operation: 

• Helicopter mode. 

• Conversion mode. 

• Cruise mods. 

Each of these modes requies a different configuration of the aircraft, i.e., a different 
range of mast angles for the rotors. Of particular interest is the conversion mode, 
since the most sig.iificant interference between the rotor wakes and the horizontal 
tailplane occurs during this mode. 

A typical flight profile for a typical tilt rotor aircraft, the NASA/Army/ 

Bell XV-15, is shown in Figure 19 of Reference 48 for the standard case of 13,000 lb 
gross weight, trimmed level flight at sea level, and aft center of gravity. Figures 14, 
15, 16, and 17 show the variation of rotor speed, rotor attitude, rotor thrust coeffi- 
cient, end blade collective pitch with forward speed for this standard case as given in 
the reference. Note that rotor incidence is related to rotor attitude by the relation 



-64 - 


ABROBPACB BVBTBMB, INC. • onb viNB anooK park • purlinoton, MASBACHUsarra oieoa • (ai7) 078- 


ROTOR ATTITUDE Op. DEGREES ROTOR RPM 



Figure 14. Rofor rpm versus Airspieed. 



Figure 15. Rotor Attitude versus Airspeed. 
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THRUST COEFFICIENT 




Figure 17. Rotor Collective Pitch versus Airspeed. 
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It is desired to estimate the values of colleci-ive pitch and rotor incidence for 
other flight cases of par+'cular interest in the study of aerodynamic interference effects. 
This can be done on an incremental basis with respect to the standard case. 

Consider the effect of varying gross weight, altitude, and climb angle on 
rotor attitude and collective pitch: 

• Gross Weight 

During conversion in the standard flight case, it is assumed that the aircraft 
is operated near its maximum lift-drag ratio and with a lift coefficient near unity. These 
conditions imply a drag polar of the form 

Cp = a(1+C2) (1,7) 


Now consider the change of drag coefficient resulting from the change of lift required by a 
change in gross weight: 

Cp + ACp = o[i+(Cl+AC^)2] 

-- a 1^1 + (Cj +2 ClAC^ 

« a(l + C^) + 2aC|^ AC|^ 


AC, 


2aC^^ACL 


(119) 




( 120 ) 
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It is seen that for original lift coefficients close to unity, the percent change in drag 
coefficient is nearly equal to the percent change in lift coefficient. Therefore, to 
maintain equilibrium along the flight path, only a change in thrust, i.e., collective 
pitch, is required. No change in rotor attitude is necessary. 

The required change in thrust coefficient is 


AC^ = 


\'^i 


( 121 ) 


• Altitude 

Consider the change of drag coefficient resulting from the change of lift 
coefficient required by a change in altitude, i.e., air density. As for changes in 
gross weight, the percent change in drag is nearly equal to the percent change in lift. 
Therefore, only a change in thrust, i.e., collective pitch, is required; 


A C 


T 



( 122 ) 


• Climb Angle 

During climb or descent at angle y to the horizontal, there is an incre- 
mental drag force given by 


AD = W sin Y (123) 

where 

W = aircraft gross weight 
Y = angle of climb (+) or descent (-) 
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Since there Is only a negligible change in lift required during clirr.b, a change in rotor 
attitude is required to maintain equilibrium along the flight path. In addition, a change 
in thrust, i.e., collective pitch, is required. 


Consider the requirements for equilibrium along the flight path. The incre- 
mental component of thrust coefficient along the flight path is given by 

A = (Cj + A Cy) cos (ofj^ “ ^ 

The corresponding incremental component of thrust coefficient perpendicular to the flight 
path is given by 


A C 2 = (Cj + A Cy) sin (of|^ - Acti^) - Cj sin 


( 125 ) 


Since the change in is negligible, set A - 0: 


.*. + A Cy - 


Cy sin 

sin (oTp - A Op) 


( 126 ) 


AC, = C, 


Sin O', 


sin (ap - Aa ) 


( 127 ) 


Substituting into the expression for 


Sin O', 


^T 


tan (af|^ - Aoj^) 


- Cy cos 


( 128 ) 


.*. tan (Op - Aop) = 


sin Or 


L .. 


+ COS Or 
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where 


AC^ = 


VV sin Y 
PttR^ 


(130) 


Equalion f129) can be solved for - Aoj^) and then Equation (126) yields AC-j-. 

The change in collective pitch A9 q required by a change in thrust coeffi- 
cient ACy can be obtained from the thrust expression (Eauation (108)). 


^ _ U U 


OQ 

2 


sin 0 


0.75 


1 . 3 2\ 

1 + - p 1 

2 / 


1 


+ - p 9 ^ (cos 9rt Tc + \ sin 9 
2 s 


0.75 


0.75^ 


— cos 9 

2 


0,75 


(131) 


Setting Cj - Cy + ACy and 9 q y. = Qq y^ + A9q y^, and assuming small A9 q y^, 
the incremental thrusi coefficient is 


A. /- <7 0 

ACt » — 
' 2 


cos 0 


0.75 


14-3 2 

1 + - p 


sin 0 


0.75 , 1 . Q J a . ^ Q 

^® 0.75 ' ‘=°^® 0.75 


(132) 


and the required ’■elation between ACy and A9 q y^ is obtained. 

Equations (121), (122), H27), and (129) were applied to the calculotion of 
the required thrust coefficients nd rotor attitudes in the following equilibrium flignt 
conditions; 
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• W - 10,000, 13,000 and 15,000 lb, y = C^, h = 0 ft 

• W = 13,000 lb, Y - 0®, h - 10,000 ft 

• W = 13,000 lb, Y = ±10°, h = 0 

Tbe results are shown in Figure 15 and 16. 

Equation (13?) was used to calculate the required values of collective pitch 
for the various flight conditions. Results v/ere shown in Figure 17 and in Table 5 below: 


Table 5. Trim Values of Rotor Attitude and Collective Pitch. 
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SECTION 9 


RESULTS 

The formulation presented in Sections 6 and 7 was incorporated into a 
computer program to investigate the aerodynomic effects of tilt rotor aircraft. In partic- 
ulor, the problem ot the rotor wake inte'-action effects for the XV-15 aircraft at various 
flight conditions was studied. Typical numerical results, showing the feasibility of the 
present method, arn discussed in this section. 

The geometry of the ti|r otor aircraft used is shown in Figure 18. Due 
to the preliminary nature of this study, the aerodynamic interference effects on the 
horizontal tail are assume J *o be caused principolly by the rotor wake, v/ith nacelle 
aerodynamics con’^ributions being cf a secondary nature. Therefore, the geometry of 
the engine nacelles has been omitted as seen in Figure 18. Moreover, an infinite- 
blad’ rotor, i.e., cctua'.or disk, is used. 

I'o St jay the convergence of the numerical method, the case of the tilt- 
rotor aircraft without ^e aerodynamic interference effects on ;K . tail was Investigated. 
Specifically, .. wi:.g-tail configuration and an actuator disk were treated separately and 
the-’, combined for anaivsis at the hover condition. Also, the problem with interference 
effects was demonstrated for one flight condition in the conversion mode. 

9.1 CASE W!"H NC INTERFERENCE 

Convergence studies v.ere conducted to identify ..ie approximate number of 
aerodynamic elements required. In addition, the overall vehicle lift and pitching 
moment coefficients were obtained. 

9.1.1 CONVERGENCF 

The results af a convergence study on a w.ng-tail configuration at an angle 
ot p*^tcck of 4.5° and airspeed of 80 knots without rotors are given in Figures 19, 20, 
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ACTUATOR DISK 


Figure 18. Geometry of a Till Rotor Aircraft. 
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and 21 . Since only longifudinal siability was considered, aircraft symmetry between 
the left and right sides was used. Hence, the number of elements chosen for the right- 
hand side of the oircra^t was: 


Description 

Number of Elements 

Ron 

Symbol 

Wing 

Tail 

Total 

1 

A 

32 

32 

64 

2 

o 

72 

72 

144 

3 

0 

98 

98 

196 


The variation of the potential difference A cp on the wing and tail sections of the tilting 
proprotor research aircraft at the root chord section (2''/b = 0) and the trailing edge 
(x/c = 1) is presented in Figure 19. In oddition, the corresponding upper and lower 
pressure distributions as well as the lift distributions along the wing and tail root chords 
are shown in Figures 20 and 21, respectively. Because of the effect of the wing wake on 
the horizontal tail when the wing ond tail are nearly coplanar, the pressure and the lift 
on the tail converged slowly. For this wing-tail configuration, 196 finite elements 
were determined to be sufficient for convergence. 

A convergence study was also performed for a single actuator disk. Results 
for 70, 96, and 126 elements indicated that convergence was achieved for the actuator 
disk model with 70 elements. To avoid possible inaccuracies in the total tilt rotor air- 
craft problem, 96 elements were chosen for the actuator disk representation. 

Results for the total vehicle problem, i.e., the wing-tail configuration with 
actuator disk models, are presented in Figures 22, 23, and 24. For the purpose of this 
sfudy, the fuselage was assumed to contribute insignificantly to the aircraft lift and 
pitching moment and, therefore, it can be neglected. Similarly, the contribution of the 
vertical tail to the l''ngitudinal dynamics can also be neglected. To preclude Interference 
effects on the horizontal tall, the case of the aircraft at hover condition, that is, with 
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the nacelie incidence at 90°, was considered. The number of aerodynamic elements 
used in the study is summcrized as follows: 


Description 

Number of Elements 

Run 

Symbol 

Wing 

Tail 

Actuator Disk 

Total 

1 

A 

32 

32 

32 

96 

2 

o 

50 

50 

96 

196 

3 

0 

70 

70 

96 

236 


Tne potential difference, A cp, and pressures on the upper and lower surfaces of the wing 
and tail are given in Figures 22 and 23. Again, the results show that the convergence 
of the pressure distribution on the tail section is slow and this is due to the presence of 
the wing wake on the tail. In addition, the co'^responding lift distribution on the wing 
and tail is presented in Figure 24. 

9.1.2 LIFT AND PITCHING MOMENT 

Comparison of the lift distributions in Figures 21 and 24 shows that the 
presence of the rotor produces an increase in the wing lift and a decrease in the tail 
lift distribution. The wing-tail configuration without proprotor at cn angle of attack of 
4.5° and airspeed of 80 knots was found to yield very little lift as shown in Figure 25(a). 
Also, a negative pitching moment (nose down) was obtained for this configuration 
(Figure 25(b)). In evaluating the pitching moment about the oircraft center-of-gravity 
position, some basic assumptions were used wh’ch affect the accuracy of the results. 

The thrust vector was assumed to act through the center of gravity. In addition, only the 
lift and induced drag on the wing and tail were used to determine the value of the pitch- 
ing moment; that is, viscous drag (drag due to skin friction) was not included in the 
analysis. This value can be obtained from wind tunnel tests of a tilting proprotor re- 
search aircraft, ond must be included if true trim conditions are to be calculated. 


-76 - 


ABROBPACB ■YCTBIWB, INC. • one vine anoOK pamk • sunuNOTOw, MASBACHuaETTS ciaoa 


7 ) a^a- 



The lift coefficient for on actuator disk in hover with angular speed of 
600 rpm is shown in Figure 26. The case of an actuator disk with no central core 
vortex, i.e., 6cpg = 0, is indicated by the square symbols. This represents the presence 
of the vorticity shed by the rotor blades near the tip of each blade, that is, the rim 
wake of an actuator disk. As expected from discussion in Section 6, the results show 
that no lift was obtained. In order to produce lift on the actuator disk, the formulation 
must include a central vortex, i.e., A?Pg / 0, and this is treated in Subsection 6.1 . For 
this case, the vorticity shed by each blade near the blade root as we* I as the tip was 
considered. To be exact, the wake is composed of infinite vortices emanating .rom the 
tip of the blade and one additional vortex emanating from the center of the disk. The 
results for this case are denoted by circles in Figure 26. Furthermore, no pitching 
moments about the actuator disk were obtained for the hover case. 

Next, the problem of the wing-tail-rotor configuration with nacelle incidence 
at 90° was treated. The total lift and pitdiing moment around the aircraft center o; 
gravity are presented in Figure 27. For this flight configuration, the presence of the 
additional central vortex on the actuator disk yielded greater lift and larger nose-down 
pitching moment on the tilting proprotor aircraft. Again, it it emphasized that the 
viscous drag of the aircraft was not used in the computation. 

9.2 INTERrERENCE EFFECTS 

A study of the interference effects was conducted for the tilt rotor oircraft 
In the conversion mode of flight with nacelle incidence set or 60°, airspeed of 100 
knots, rotor angular speed of 565 rpm, and a tail incidence of -1 .5°. For this 
flight condition, the horizontal h-il was imbedded in the rotor wake. The variation 
with angle of attack of the total vertical forces on the aircraft is given in Figure 
28(a). Vertical eq"ilibrium is achieved when the sum of the vertical forces is zero. 
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For this configuration, vertical equilibrium was obtained with an angle of attack setting 
at 2.3°. In addition, the analysis of this same problem without the implementation of 
the additional central vortex on the actuator disk yielded the required angle of attack 
setting to be about 2°. 

The effect of rotor wake interference on the tilting propiotor aircraft at this 
flight configuration is presented in Figures 28(b) and (c). ("hese results are for the 
aircraft in the untrimmed condition and represent the change in the lift and pitching 
moments of the aircraft with interference effects relative to the aircraft with no rotor 
woke interference effects. As seen in Figure 28(h), the total aircraft lift coefficient 
was significantly decreased by the presence of interference effects. The reduction in 
lift was less for the actuator disk with a central core vortex (Cj = .0064). The lotal 
oircrafi pitching moment in the nose down direction was significantly increased by the 
interference effects (Figure 28(c)). As seen in the figure, the addition of the central 
vortex on the actuator disk yielded a small difference in the overall pitching moment 
for this flight condition. 
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(b) Wing, Spa?, wise Direction, Trailing Edge 


Figure 19. Potential Difference, Acp, Wing-Tail Combination, 

No Rotor, Airspeed = 80 kr, a = 4.5°. 
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(c) Tail, CHordv/ise OirecMon, Root Chord 
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(d) Tail, Spanwlse Direction, Trailing Edge 

Figure 19, Potential Difference, Acp, Wing-TaU Combination, 

No Rotor, Airspeed =■■ 80 kt, a = 4.5°. (Concluded) 
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Figure 20, Pressure Distribution, Wing-Tail Combination, 
No Rotor, Airspeed = 80 kt, a - 4.5°, 
Chordwise Di'ection, Root Chore . 
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(c) Tail, Upper Surface 



Figure 20. Pressure Distribution, Wing-Tail Combination, 
No Rotor, Airspeed = 80 kt, a - 4.5°, 
Chordwise Direction, Root Chord. (Concluded) 
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(a) Wing 



(b) Tail 


Figure 21. Lift Distribution, Wing-Tail Combination, No Rotor, 
Airspeed = 80 kt, a = 4.5°, Chordwise Direction, 
Root Chord. 
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(a) Wing, Chordwise Direction, Root Chord 



(2y/b)wiN6 

(b) Wing, Spanwise Direction, Trailing Edge 


Figure 22. Potential Difference, A Wing-Tai! -Rotor Configuration, 
Airspeed = 80 kt, a - 4.5°, Nacelle Incidence = 90°, 
Angular Spew’d = 565 rpm. 
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(d) Tail, Sprnwise Direction, Trailing Edge 

. igjre 22. Potential Difference, Acp, Wing-Ta!i-Rotor Configuration, 
Airspeed = 80 kt, a - 4.5°, Nacelle Incidence = 90®, 
Angular Speed = 565 rpm. (Concluded) 
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figure 23. Pressure Distribution, W'rg-Tail-Rotor Configuration Airspeed = 80 kt, 
01 = 4 , 50 ^ Nacelle incidence = 90°, Angular Speed - 565 rpm, 
Chordwise Direction, Root Chord. 
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(d) Tail, Lower Surface 


Figure 23. Pressure Distribution, Wing -Tail -Rotor Configuration, Airspeed = 80 kt, 
a = 4.5°, Nacelle Incidence = 90®, Angular Soeed = 565 rpm, 
Chordwise Direction, Root Chord. (Conclude 
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NUMBER OF ELEMENTS 
(a) Lift G>efficient 



(b) f-itching Moment Coefficient 


Figure 25. Aerodynamic Coefficients as Function of Number of 
Elements, Wing~Tail Combination, No Rotor, 
Airspeed = 80 kt, a = 4.5°. 
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Figure ?,6. Lift Coefficient on Actuator Did< as Function of 
Number of Elements, Angular Speed = 600 rpm. 
Mast Angle = 90®. 
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(a) Lift Coefficient 
NUMBER OF ELEMENTS 



(b) Pitching Moment Coefficient 

Figure 27. Aerodynamic Coefficients as Function of Number of Elements, 
Wing-Tail-Rotor Configuration, Airspeed = 80 kt, a =4,5°, 
Nacelle Incidence = 90°, Angular Speed = 565 rpm. 
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Figure 28. Interference Effects for Conversion Mode, Wing-Tail -Rotor 
Configuratioii, Airspeed = 100 kt. Toil Incidence = -1.5°, 
Nacelle Incidence = 60°, Angular Speed = 565 rpm. 
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Figure 28. Interference Effects for Conversion Mode, Wing-Tail-Rotor 
Configuration, Airspeed = 100 kt. Tail Incidence = -1 .5°, 
Nacelle Incidence = 60°, Angular Speed = 565 rpm. 
(Concluded) 
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SECTION 10 

CONCLUSIONS AND RECOMMENDATIONS 

The theoretical methodology and corresponding numerical procedure for the 
evaluation of rotor -aircraft aerodynamics have been formulated. Results were obtained 
for a wing-body-rotor configuration of the XV-15 tilting proprotor aircraft at various 
flight conditions and they are presented in Section 9. Concluding remarks on the present 
study are summarized in Subsection 10.1 . Based on the results obtained, the areas 
recommended for future resecrch are identified in Subsection 10.2. 


10.1 


CONCLUSIONS 


The purpose of this investigation was to demonstrate the feasibility and flexi- 
bility of a computer method to study the rotor-aircraft aerodynamics with special atten- 
tion given to the rotor wake interference effects on tilting proprotor research aircraft. 

Due to the preliminary nature of the study, the geometry of the fuselage, engine nacelles, 
and vertical stabilizers were not included in the analysis. Furthermore, the presence 
of the viscous drag on the aircraft was not included in the pitching moment calcula- 
tion. 

The numerical technique used in this investigation indicated a good rate of 
convergence with the possible exception of the pressure on the tail. This was Oue to the 
effect of the wing wake on the tail. The subject of coplanar wake interference is 
treated in great detail in Reference 49. For the flight configurations explored in this 
study, the effects of the rotor wake interference on the XV-15 tilt rotor aircraft yielded 
a reduction in the total lift and on increase in the nose-down pitching moment. 

In order to extend the existing theoretical formulation to include actuator 
disk aerodynamics, new concepts were developed. As shown in Appendix B, the tilt 
rotor aircraft problem yielded a system of algebraic equations with zero determinant. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


95 


mamma tumi waa nma 


ABROBf»ACB ■VBTBMa, INC. • ONB VINB anOOK PARK • aUPLINOTON, MASBACHUaBTra 01803 • (817) 978-7017 


..,1 



It must be emphasized that the problem of zero determinant is encountered whenever the 
wake separates the flow field into two regimes and is not peculiar to the actuator disk. 
The problem also occurs with the aerodynamics of blunt fuselage configurations as well 
as with aerodynamics of buildings. Also, in solving <’he problem with the actuator disk, 
the strength of the branch wake, ^ 9 ^, may assume any value and cannot be determined 
explicitly in general. In Appendix C it is shown that A cpg can be prescribed as a func- 
tion of rotor thrust. In addition to exploration of the problem being singular, new 
boundary conditions for the actuator disk were formulated and discussed in Section 6 . 

10.2 RECOMMENDATIONS FOR FUTURE RESEARCH 

The feasibility and flexibility of a computer method to study rotor wake inter- 
ference effects on the XV-15 aircraft have been demonstrated. The computer program 
can easily be extended for analysis of more complex problems such as helicopter and 
other VTOL configurations. It is emphasized that the only required inputs to the program 
are the locations of the corner points of the aerodynamic elements and certain specifica- 
tions to the problem, i.e., Mach number, rotor angular speed, angle of attack, and 
rotor thrust for the actuator disk. In addition, the geometry of the wake is automatically 
generated. It is also of particular importance to the user that the program is extremely 
easy to use because familiarity with the specific method is thus not required. Present 
outputs available from the computer program include the aerodynamic coefficients, i.e., 
lift, induced drag, and pitching moment as well as the pressure and velocity flow field 
around the configuration. 

It is emphasized that the use of an actuator disk model for the aircraft rotor 
allows the study of aerodynamic interference effects on the XV-15 tilting proprotor 
aircraft, but does not allow the prediction of the lift on the rotor. The Green's formu- 
lation for the actuator disk requires the strength of the branch wake on the disk to be 
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prescribed as an Input and this strength is proportional to the lift or thrust on the rotor. 
Hence, to obtain the lifr on the rotor, a formulation with rotor blades can be developed 
and included in the analysis. 




More detailed study on the XV-15 tilt rotor aircraft can easily be undertaken 
to obtain and assess the following: 

• Effects of a rotor-blade model on the aircraft aerodynamic interference. 

• Effects of vehicle configuration and speed on the longitudinal 
trim. 

• Sensitivity of the aerodynamic interference effects to variations in 
aircraft gross weight and power. 

• Effects of aerodynamic interference on the aircraft stability and 
handling qualities and characteristics. 

Other areas for future investigation include: 

• The case of fully unsteady incompressible flow can be analyzed. This 
requires a step-oy-step analysis in time: at each step a new geometry 
of the wake is obtained from the tangency condition, while the shed- 
vorticity intensity is known from preceding time histories. Also, the 
actual geometry of the tilt rotor aircraft can be implemented. Since 
the problem is unsteady, the actuator disk model is not required. 

• The effect of the compressibility can also be included with special 
attention given to the analysis of the singular blade problem. This 
reaiire'. an extensive modification of the procedure for incompressible 
flow. 
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APPENDIX A 

AERODYNAMIC POTENTIAL PROGRAM SUSSA ACTS 

Computer program SUSSA ACTS (Steady and Unsteady Subsonic and Super- 
sonic Aerodynamics for Aerospace Complex Transportation Systems) implements the 
method of Morino (References 27 to 29 and 43 to 46). This program is very flexible, 
simple to use, and modular. The geometry is defined in one section, while a second 
part evaluates the normal wash, the porentiai, the jressure coefficient, and the aero- 
dynamic coefficients (lift, moment, and generalized-force coefficients). Moreover, the 
user need not be familiar with the aerodynamic portion of the program, unlike other 
sop'^isticated aerodynamic programs in which the choice for the combination of various 
panels (sources, doublets, vortices) is an art which requires considerable understanding 
of the method. Another advantage of SUSSA ACTS is that the paneling used for the 
aerodynamics is completely arbitrary and, therefore, may coincide with the one used 
for the structural analysis. These attributes of the program make the method extremely 
efficient from a practical paint of view. Finally, the method is both accurate and 
fast, despite the fact that no effort has yet been made to minimize the computation time 
(see References 29, 44, and 46). Other important features of SUSSA ACTS are: 

• Inclusion of control surface mode shapes (i.e., surface rotation). 

• Option to compute aerodynamic forces caused by turbulence. 

• Computation of downwash from arbitrary mcdes. 

c Geomelry preprocessor for conventional wing-body-tail 
combination , 

Option to compute low-frequency (quasi -steady aerodynamic 
coefficients). 

A flow diagram of SUSSA ACTS is presented in Figure A-1 . The sub- 
routines which perform the indicated functions are denoted in the appropriate blocks. 
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Figure A-1 . Program SUSSA ACTS. 
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Basically, SUSSA ACTS performs \yto separate functions: one is the definition of the 
aerodynamic paneling on the aircraft, and the second is the computation of the aero- 
dynamic potential and pressure distributions and the determinoiion of the aerodynamic 
force coefficients. The aerodynamic paneling can be accomplished manually by inputing 

the coordinates of the corner points of all selected aerodynamic panels, or by using a 

I 

geometry preprocessor. A geometry preprocessor for conventional wing-body-tall con- 
figurations is incorporated in SUSSA ACTS (Reference 47). Given the aircraft con- 
figuration geometry, the preprocessor computes the corner point locutions of required 
aerodynamic panels for the wing-body-tail configuration. Required inputs include: 

• Fuselage length. 

• Fuselage diameter (assumed circular) - center section. 

• Fuselage shape factors for nose and tail sections. 

• Wing span. 

e Wing sweep angle. 

• Horizontal and vertical tail ^ans. 

• Horizontal and vertical tail sweep angles. 

This routine eliminates the requirement for Hie selection of aerodynamic paneling by 
the user. 

The second section of SUSSA ACTS implements the aerodynamic potential 
computation method of Reference 44. These computations require the aerodynamic 
panel corner point locations and input rigid-body and flexible-body mode shapes, 
control surface rotations, and turbulence. The aerodynamic potential is computed 
first, the pressure coefficient distributions are then determined, and the aerodynamic 
force coefficients are computed. 
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APPENDIX B 

INTEGRAL EQUATION FOR ACTUATOR DISK 


A general Integral equation for actuator disks is presented in Section 6. In 
this appendix a few mathematical details which are important for the solution of the 
Integra equation are discussed thoroughly. For simplicity the emphasis is placed on 
o . a; i /mmetric problem, although all the main conclusions are valid for a general 
'configuration as well. 

B. 1 INTEGRAL REPRESENTATION FOR ACTUATOR DISK 


The integral representation of the potential for actuator disks is presented in 
Sectio 6. Noting that since the potential difference A cp = A cpg is constant on the 
i>ranch wake. Equation (69) may be rewritten as 


4tt 


^ 1 

9n \ r / 3n 


Ed " ' 


dE 


(B-1) 


^D “B 

If P* is Eq, the function E(P) assumes the value E* = 1/2 and 
Equatiori (B-1) reduces to an integral equation. As mentioned in Section 6, in solving 
this integral equation, Acpg may assume any value. Therefore, in the actuator disk 
formulalion, A cpg (which is approximately proportional to ri's thrust) is nol an output 
but an input to the problem. This indicates that the actuator disk integral equation is 
singular, since solution is not unique. This point is analyzed further in Subsection 
B.3. 

in addition, even if A cpg is treated as a constant, the integral equation is 
s-.il singular. The problem exists any time when the wake separates the flow field into 


B - 1 
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two separate regions, one inside and one outside the wdce. It should be noted that 
this is the only reason for this problem, and therefore is not peculiar to the actuator 
disk. This problem occurs in the formulation for bodies with a blunt back such as a 
fuselage with truncated base. A similar prc^lem occurs in the formulation for building 
aerodynamics. This point is analyzed in detail in Section B.2. In order to emphasize 
that the problem is not peculiar to the actuator disk formulation, A cp is assumed to be 
equal to zero (i.e. , there is no vortex emanating from the root of the blades), and 
this is presented in Section B.2. The extension to the case with Atpg different from 
zero is considered in Section B.3. 


B.2 ACTUATOR DISK WITH = 0 

B 

In order to show that the integral equation for the actuator disk is singular 
even if A 9 g is prescribed, it is convenient to obtain Equation (B-1) a new way. Assume 
Acp = 0. Consider Figure B-l(c), and note that the value E = 1 outside the actuator 
disk may be obtained as the sum of the function E. for inner flow (i.e., inside the 
wake) and the function E^ for the outer flow (i.e., outside the wake) or 


E(P) = E.(P) + E^(P) 


(B-2) 


yields 


and 


Applying Green's theorem to the regions shown in Figures B-l(a) and B-l(b) 




dZ. 




dE. 


(B-3) 


(B-4) 
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Note that 


n. = n ; n = n on Er. 
I o D 


n. = -n; = n on Ej^ 


(B-5) 


Thus, adding Equations (B-3) and (B-4) and using Equations (B-2) and (3-5) yields 

4ttE(P*)cp(P*) = f-j "‘P ^ 

.JJa,a(i) d. 


where Acp = cp^ " 9j • Note that Equation (B-6) coincides with Equation (B-1) with 
A cpg =0, In order to show that the operator in Equation (B-6) is singular, it is suffi- 
cient to show that the homogeneous equation (i.e,, with 9cp/S n = 0) 


4ttE(P*)9(P*) = + J/ 


dE (B-7) 


has a nontrivial solution. In order to show this, consider the function 


9 (P) = 0 P outside E. 
NT ' 

= 1 P inside E. 


If 99/9n - 0 (homogeneous problem), this function is a trival solution 
V .e., 0 = 0) of Equation (B-4), whereas for Equation (B-3) one obtains (note that 

Ei«PNT ” ^i> 


E. • ' 

I 
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which Is identically satisfied, since (indicating with Q the solid angle) 


ijdZ. sCi = 0 outside S. 

= 2tt on Z. 

I 

= 4n inside Z. (B-10) 

Since cpj^y is a nontrivial solution for both Equations (B-3) and (B-4) with 3cp/9n = 0, 
it is a nontrivial solution for their sum, Equotlon (B-7). Hence, the operator in 
Equation (B-1) (with A cpg =0) is singular; therefore, if cp satisfies Equation (B-1), 
then 

cp = 9 + Ccpj^^ (B-11) 

(where C is an arbitrary constant) also satisfies Equation (B-1). 

Finally note that in the numerical formulation. Equation (B-1) is replaced by 
Equation (91). However, the doublet Integral represents solid angles and the total 
solid angle is evaluated exactly with the use <• f rh- hyperboloidal elements. Thus the 
discrete form of Equation (B-10) Is still valid e..actly. Therefore, even the discrete 
system is singular, i.e., the determinant of the system given by Equation (91) is equal 
to zero. This implies that the vector 

{»r"} = {'pntC’i'} 

is a nontrivial solution for Equation (91) with right-hand-side equal to zero, or that 
if is a solution to Equation (91), then 

{t.) = {$,} + C { 9 ^^} (B-13) 

(where C is on arbitrary constant) is also a solution to Equation (91) (i.e., if the 
solution exists and is not unique). 
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B.3 ACTUATOR DISK WITH Acpg / 0 

The operator in Equation (B-1) has an additional nontrivial solution; however, 
this solution does not satisfy Equation (B-6) which is obtained from Equation (B-1) by 
dropping the branch wake. The implications of this feature are discussed later in the 
section. An explicit solution for Ec/jotion (B-1) is given for an actuator disk in axial 
flow, i.e., for an oxi mmetric case. Then it is shown that the general case has a 
solution with features similar to the ones of the axisymmetric case. 

The geometry for the axisymmetric case is given in Figure B-2, where again 
it is shown that the complete flow (Figure B-2(c)) may be thought of as the superposition 
of the flow outside the wake (Figure B-2(b)) plus the one inside the wake (Figure B-2(a)). 
The shaded area indicates the branch wake (which does not exist for the outer flow. 
Figure B-2(b)). Applying Green's theorem to the geometries of Fig^ es B-2(a), B-2(b), 
and B-2(c) yields, respectively 



+ A 



dE 


(B-1 4) 


where Acpg is the (constant) potential discontinuity on the branch-wake surface. Eg 


4rE^(P*)cp(PJ = 



(B-15) 


B-6 


ABROe'*ACH avaTaMISi INC. • on* vini aMaoK pa^k • BUPLiNaTON, MABSACHuarrra oibo3 • (ai7) a7a-7Bi7 



(c) Flow Outside the Actuator Disk 


Figure B-2. The Geometry for Axisymmetric Actuator Oidc Problem. 
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and 


4itE(P.)<|,(P.) = -(ft 

" L«n \r/ 3n\r/J 

3n 


d£ 


+ Acpg 


JJ- 

J ""B 


L] dS 


(B-16) 


Nofe that Equation (B-16) may be obtained by adding Equations (B-14) and (B-15). 

Consider the nontrivial solution 

9b(P) - ^ for P inside 2. (-n < 9 < n) 

= 0 for P outside 2. (B-17) 

where 6 is defined in Figure B-3(a). Note that, by definition of 6, the branch wake Is 
located at 9 = ±tt . This is not a limitation since, according to Equation (B-14), it is 
always possible to add a constant to (p„. Note that in any event 


A?b = 


(B-18) 


If 39/3n = 0 (homogeneous equation), cpg satisfies trivially (i.e., 0=0) 
Equation (B-15) whereas for Equation (B-14) one obtains 


4l-E(P.)<p(P.) = ^ 9 


^ ll dE, 


_ 3n. \ r 

2 . 1 ' 

I 




(B-19) 
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(b) Definition of 9^ OF POOR QUALITY 

Figure B-3. Actuator Disk Branch Wake Potential Geometry for Nontrivial Solution. 
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Note that (see Figure B-4) 


A L dE, = --L (r ..„.)dE = 

an. \rj ' ,3 ‘ ,2 


dZ 


= dCl 


(B-20) 


is the solid angle for dl as seen from P* (positive since (n -r) < 0)). In ordi ’ to 
evaluate t. e first integral it is convenient to introduce the varioble 

9| =9-9* 

The function cpg cs a function of 0^ is given in Figure B-5, Limiting 0^ to the interval 
(“TT , Tt) the function cpg is given by 


cPb(P) = 9^ + 0,^ - 2ttH(9) 


where (ossuming 9* > 0 for simplicityt) 
H(9,) = 0 
= 1 


-rr<9^ <n - 9* 


n - 9^<9j <n 


Thus, Equation (B-20) may be rewritten as 


4ttE*(P*) 9(P*) = //[9^ +9*-2nH(9)| da 
-2TT JJdng 


(B-21) 


(B-22) 


(B-23) 


where dQg is negative since n^ ' r > 0. 


^This hypothesis is not restrictive in view of the axisymmetry of the problem. 
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Figuro B-4, Solid Angle 
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Next, note that 

II e, dn = 0 (B-24) 

r. 


because of antisymmetry of the integrtnd, (see Figure B-3(b)). Then, since dflg is 
negative, 

II H( 6 )dn ^ If ~ ■ If 

I. E* E' 

I I D 

Thus, Equation (B-25) reduces to 


Since 


4 TTE*(P)cpg(PJ = e* II dn = e, 4 TTE(Pj 

Ei 


11 


dn = 4 tt 
= 2n 


= 0 


if P* is inside E. 
if P* is on E. 

I 

if P is outside E. 

I 


(B-26) 


(B-27) 


Equation (B-27) is equivalent to Equation (B-17). Therefore, cpg is a nontrivial solution 
of Equation (B-16) with 3<p/9n = 0, which indicates that the operator in Equation (B-16) 
is singular. 

Finally, the above results imply that the solution of Equation (B-6) (without 
branch wake) also satisfies Equation (B-1), since that solution has A 9 g = 0. Therefore, 
not only is the determinant equal to zero, but A cpg cannot be determined by Equation 
(B-1). In order to have a solution with A cpg / 0, the value of A cpg must be prescribed 
and the branch wake integral must be treated as a known quantity. This implies that 
the boundary conditions on 8cp/3np (Equation (6-11)) are not sufficient to determine the 

problem, but the value of A cpg (proportional to the lift) must be presc. ibed as well. 
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The above solution could be obtained by noting that the circulation is equal 
to zero for o contour outside (since the contour may be reduced to zero without 
intersecting any singular line or surface) whereas the circulation is constant and differ- 
ent from zero if the contour is inside and simply connected with the vortex line 
emanating from the center of the actuotor disk. In other words, for a nonaxisymmetric 
configuration the nontrivial solution is still equal to zero outside whereas there is 
a constant discontinuity on the branch wake. 
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APPENDIX C 


EVALUATION OF THRUST ON TtlE ROTOR 


The rotor thrust perpendicular to the tip path plane can be determined From 
momeni^m considerations. For simplicity, only o single rotor blade is explicitly 
included in the analysis, although the result is valid for multi-blade rotor configura- 
tions as wel I . 

Conservation of rriomontum states that the time rate of increase of momentum 
within a fixed coritrcl volume, V, is equal to the rate ot which momenti'm is changing 
within V pius Hie net forces acting on the fluid within V. For the case of steady 
flow w'th rt-3gligible viscous and body forces, the conservation of momentum is written 


p ndS + 


pV (V • n) dE =0 






(C-1) 


where the first term represents the forces acting on the control surface due to the normal 

pressure, p , and the second term represents the change of momentum of the fluid 

in the control volume. Note that n is the unit normal directed positively out from the 

surface E|^ + Eg. The control surface Ep defines an actuator disk which encloses the 

rotor blade. As shown in Figure C-1, Sp is the contour of the rotor blade and, therefore, 

K 

there is no momentum flux through it. 

The fluid enters the control volume at station 1 and leaves at station 2. 

Note that p^, V^', and are average values of the pressure, velocity, and density, 
respectively, upstream of the rotor. Similarly, p 2 , ^2' ?2 average values 

downstream of the rotor. The lateral side of Ej^ can be chosen along streamlines and 
momentum flows through Eq only at stations 1 and 2; hence, any flux through the 
lateral side of Eq is assumed to be negligible. &i|uation (C-1) can now be written 
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* § pnd 

"d "r 

+ ^ pV (V • n)dEQ = 0 (C-2) 

The flux of momentum, m, through the control surface is 
m = ^ pV (V • n) dEQ 

= // P [(V • n), V, + (^ • n)2 V2] dE (C-3) 


Note that the continuity equation yields, in average, 
(V • n)^ = -(V • n)2 


(C-4) 


and noting that 



(C-5) 


one obtains 


(V^ - ^2) * ' = ° 


(C-6) 


Therefore, the flux through can be written as 

m = - IJ pV,- T (V, - 


(C-7) 
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Let F denote the force on the rotor 




F = - appndEj 

Er 


(C-8) 


Using Equation (C-8) in Equation (C-2) yields 

JJ(P, -P 2 )idEo-JJ PV,.T(V, -Vj)dSp 

The thrust perpendicular to the rotor tip path plane is, by applying 
Equation (C-6), 

T =J} (Pi ~P2)dEp 


( (C-9) 


(C-10) 


This is the desired expression which relates the forces acting on the rotor to the forces 
acting on the actuator disk. 

The pressure is obtained from Bernoulli's equation 


+ £.+ L Vcp*Vcp= constant 

at p 2 

it is assumed that the radial velocity is negligible when compared with 
circumferential velocity and this yields 


(C-11) 


d(p _ 

at 


d0 a^ 
dt ae 


= n ^ 
ae 


(C-12) 
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Equation (C-11) con be written as 


9 9 / 


Pi -P2 = " p 

Combining Equations (C-10) and (C-13) and neglecting higher order terms in Equation 
(C-13), the thrust on the rotor becomr*. 

.R 2n 9(cp^ -cp2> 


(C-13) 


T = -pfl [ f 
0 0 


d0 rdr 


be 


= - pfl r Tcp^ -CP2] rd 

■^0 


(C-U) 


Note that 


0 

h-'.f-t/r-wr 


= -A 


because the branch wake surface is on the lower surface of the disk, 


Combining Equations (C-i4) and (C-15), 


(C-15) 


1 2 

T = - pnR^ A: 

2 


(C-16) 


which is an expression which relates the discontinuity on the wake branch with the 
rotor thrust. In terms of the thrust coefficient. Equation (C-16) becomes 


S = 


- pCpR^nR^ ttOR^ 
2 


(C-17) 
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where 


A 9g = strength of branch wake 

R = rotor radius 
n = rotor speed 


t 

/ 
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